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ABSTRACT 

Inferences  regarding  the  character  of  the  pointwise  or  material  flow 
behavior  for  cohesionless  soils  are  drawn  from  the  results  of  a  series  of 
tests  on  a  sand  system.  Specifically  these  inferences  are  that  normality 
to  the  failure  surface  may  be  preserved  for  projections  of  the  plastic 
strain  increment  vector  onto  planes  of  constant  hydrostatic  pressure 
(ir-planes)  and  that  the  angular  deviation  from  normality  to  the  failure 
surface  for  projections  of  the  plastic  strain  increment  vector  onto  planes 
perpendicular  to  the  ir-plane  may,  as  a  first  approximation,  be  considered 
constant.  For  this  latter  inference  the  variation  in  the  approximately 
constant  angle  appears  to  be  a  linear  function  of  the  relative  density  of 
the  system  at  failure.  The  apparatus  used  to  perform  the  experiments  was 
a  standard  direct  shear  device.  Statistical  methods  were  employed  to  es¬ 
tablish  a  nest  of  failure  envelopes  in  the  force  space  of  the  applied 
tractions  with  a  parametric  dependence  on  the  relative  density  at  failure. 
Normals  to  these  surfaces  were  systematically  compared  to  their  associated 
displacement  vectors.  The  results  of  this  analysis  led  to  the  inferences 
regarding  material  behavior. 
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INTRODUCTION 

The  history  of  the  application  of  the  mathematical  theory  of  plasticity 
to  problems  in  soil  mechanics  contains  numerous  examples  in  which  an  analytic 
form,  chosen  to  represent  the  failure  envelope  for  the  material,  has  also  been 
used  as  a  yield  surface.  The  stress-strain  relations  resulting  from  this  as¬ 
sumption  require  that  a  certain  dilatation  rate  be  in  evidence  when  the  stress 
state  reaches  this  yield  surface.  Experiments,  however,  have  consistently 
shown  smaller  rates  of  volume  change  than  those  predicted.  This  anomaly  has 
given  rise  to  considerable  theoretical  conjecture  about  the  inelastic  behavior 
of  soils,  yet  the  necessary  experiments  to  substantiate  these  hypotheses  have 
been  few.  The  purpose  of  this  paper  is  to  report  the  results  of  a  series  of 
tests  designed  to  ascertain  the  nature  of  the  flow  behavior  of  sands  at  fail¬ 
ure.  Specifically  it  was  desired  to  determine  if  any  relationship  existed 
between  the  direction  of  a  normal  to  the  failure  envelope  and  the  direction  of 
the  observed  plastic  flow  rate  vector  associated  with  that  normal. 

In  accomplishing  this  objective  a  somewhat  irregular  experimental  procedure 
was  used.  The  basis  for  this  technique  is  more  fully  presented  in  the  next 
section,  but  essentially  it  involves  studying  the  phenomenon  under  investigation 
at  a  lower  level  of  abstraction  than  the  customary  stress-strain  level.  Such  a 
procedure  requires  that  only  the  force  and  the  corresponding  displacement  histor¬ 
ies  acting  on  the  boundary  of  a  soil  system  during  a  test  be  monitored  in  order  to 
extract  useful  information  about  the  flow  behavior  of  sand  at  failure.  While 
such  a  technique  circumvents  the  difficulties  associated  with  producing  and 
maintaining  a  homogeneous  state  of  stress  and  strain  in  a  soil  specimen,  it  has 
the  obvious  drawback  that  the  results  obtained  must  be  regarded  only  as  infer¬ 
ences  to  the  actual  behavior  on  the  stress-strain  level.  Although  this  draw¬ 
back  is  of  serious  consequence  to  the  formulation  of  quantitatively  predictive 
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constitutive  relationships  for  the  material,  the  paucity  of  experimental  in¬ 
formation  in  this  area  and  the  still  questionable  validity  as  to  the  homogenuity 
of  stresses  and  strains  near  failure  in  specimens  tested  in  more  sophisticated 
devices  indicate  that  the  results  presented  here  may  be  of  interest. 

The  testing  apparatus  used  in  this  investigation  was  a  standard  direct 
shear  device.  The  soil  specimen  in  this  device  is  subjected  to  certain  applied 
tractions  on  its  boundaries  and  both  these  tractions  and  their  corresponding 
displacements  can  be  measured.  The  stress  and  strain  states  within  the  speci¬ 
men  are,  however,  complex  and  unknown.  As  the  testing  procedure  for  this  appara¬ 
tus  is  comparatively  simple  and  the  speed  at  which  tests  can  be  performed  com¬ 
paratively  rapid,  the  use  of  the  device  offers  a  somewhat  expedient  method  for 
obtaining  a  large  quantity  of  data.  Increasing  the  amount  of  data  enhances  the 
reliability  of  the  inferences  drawn  from  the  analysis  of  the  results  by  the  use 
of  statistical  methods.  While  such  reasoning  cannot  in  itself  justify  the 
choice  of  the  direct  shear  device,  the  well  established  observation  that  the 
force-deformation  curves  obtained  from  direct  shear  testing  do  qualitatively 
imitate  their  triaxial  test  counterparts  offers  some  additional  justification 
for  its  adoption  in  this  testing  program. 

Conceptually  the  experimental  program  involved  the  establishment  of  a 
failure  envelope  in  a  hyperspace  of  the  applied  tractions  and  the  configuration 
parameters  of  the  system  which  best  fitted  the  observed  data.  Equally,  this 
single  failure  envelope  may  be  represented  as  a  nest  of  envelopes  in  a  space 
of  the  applied  tractions.  Normals  to  this  nest  of  failure  envelopes  were  then 
systematically  compared  to  their  corresponding  plastic  deformation  increment 
vectors  to  ascertain  if  any  relationship  existed  between  them  at  failure. 

The  subsequent  section  presents  the  historical  background  to  this  subject 
and  the  basis  underlying  the  adoption  of  the  particular  investigative  approach 
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used  in  this  study.  Succeeding  sections  contain  a  description  of  the  experi¬ 
mental  procedure,  a  discussion  of  the  testing  program,  the  presentation  of 
the  test  results,  and  the  analysis  made  of  these  results.  Lastly  the  infer¬ 
ences  obtained  from  the  analysis  are  summarized  and  discussed.  The  definitions 
of  all  symbols  used  in  the  paper  may  be  found  in  Appendix  II.  -  Notation. 

Each  symbol  will  be  explained  in  the  text  when  first  used. 

HISTORICAL  BACKGROUND 

The  application  of  the  mathematical  theory  of  plasticity  to  the  problems 

of  soil  mechanics  was  initially  made  in  a  series  of  papers  by  Drucker  and 
2 

Prager  (6)  ,  Shield  (28),  Drucker  (7),  and  Shield  (29,  30).  In  these  investi¬ 
gations  generalizations  of  the  Coulomb  failure  criterion  were  proposed  as  pos¬ 
sible  yield  surfaces  for  the  material.  As  pointed  out  by  these  writers  and 
later  commented  on  by  Drucker,  Gibson  and  Henkel  (9)  the  use  of  such  general¬ 
izations  as  possible  yield  surfaces  places  certain  restrictions  on  the  direction 
of  the  plastic  strain  increment  vector  at  failure.  One  of  these  restrictions 
requires  a  certain  dilatation  rate  to  accompany  a  stress  state  on  the  yield 
surface.  Hansen  (13)  observed  that  the  dilatation  rates  predicted  by  plasti¬ 
city  theory  were  considerably  in  excess  of  those  found  in  practice.  This  dis¬ 
crepancy  between  prediction  and  observation  has  been  substantiated  numerous 
times  in  the  literature. 

Subsequent  theoretical  work  by  Drucker  et  al  (9),  Jenike  and  Shield  (18), 
Drucker  (10,  12),  Rowe  (25,  26),  Brown  (4),  Horne  (17),  Palmer  (21)  and  this 
author  and  Paslay  (31)  has  been  directed  towards  the  establishment  of  an  inelastic 
theory  which  more  accurately  predicts  the  behavior  of  granular  media  and  thus 
accomodates  the  aforementioned  discrepancy.  One  of  the  purposes  of  performing 
the  experiments  reported  here  was  to  provide  data  which  might  be  useful  in 

2  I  .... 

Numbers  in  parentheses  refer  to  corresponding  items  in  Appendix  I.  -  References. 
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testing  the  validity  of  these  more  advanced  theories. 

In  conjunction  with  the  analytical  work  mentioned  above,  considerable 
advancement  has  been  made  in  the  field  of  mechanical  testing  of  soils.  The 
modifications  to  triaxial  devices  and  procedures  by  Rowe  and  Barden  (27) t 
Bishop  and  Green  (2),  Barden  and  Khayatt  (1)  and  others  represent  considerable 
improvements  towards  the  achievement  of  homogeneous  states  of  stress  and  strain 
within  the  triaxial  sample.  Apparati  devised  by  Kirkpatrick  (19),  Haythomthwaite 
(14),  Wu,  Loh,  and  Malvern  (33)  and  Ko  and  Scott  (20)  have  been  notable  advance¬ 
ments  in  the  technique  of  applying  a  truly  three  dimensional  stress  field.  Out 
of  this  work  and  the  contributions  of  Roscoe,  Schofield  and  Wroth  (23),  Roscoe, 
Schofield,  and  Thurairajah  (24)  and  Comforth  (5)  a  more  accurate  picture  of 
the  shape  of  the  failure  envelope  for  granular  media  and  its  dependence  on 
material  density  has  emerged.  Still,  however,  certain  doubts  persist  as  to 
the  statical  determinateness  and  homogeneity  of  the  stress  and  strain  fields 
within  the  specimens  tested  by  these  more  sophisticated  devices.  See,  for  ex¬ 
ample,  Haythomthwaite  (15,  16)  and  Broms  and  Jamal  (3).  Unfortunately  the 
greatest  uncertainty  exists  at  or  near  the  failure  of  the  specimen.  How  well 
the  various  testing  devices  are  predicting  pointwise  material  properties  rather 
than  system  behavior  in  these  limiting  regions  is  still  a  question  of  some  valid¬ 
ity.  As  it  was  precisely  in  this  neighborhood  of  failure  that  information  con¬ 
cerning  the  flow  behavior  was  sought,  the  experimental  procedure  devised  for 
this  investigation  attempted  to  evade  the  issue  of  stresses  and  strains  by 
studying  the  phenomenon  at  a  lower  level  of  abstraction. 

Such  a  technique  rests  on  a  certain  concept  developed  within  the  field 
of  plasticity  by  Prager  (22)  and  furthered  by  Haythomthwaite  (14)  and  Drucker 
(11).  Consider  a  material  having  a  well  defined  yield  surface  in  stress 
space.  A  two  dimensional  projection  of  this  surface  is  shown  as  the  convex 
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closed  curve  f  (o^j)  =  K  in  stress  space  of  Fig.  1.  If  a  homogeneously 

stressed  sample  of  the  material  is  tested  such  that  the  loading  point  moves  from 
some  point  interior  to  the  yield  surface  along  any  path  to  some  point  on  the 

yield  surface  (path  OY,  Fig.  1)  then  plastic  deformation  ensues  on  reaching 

Y  P 

.  Superposing  the  axes  of  a  plastic  strain  space  on  those  of 

the  stress  space  ,  the  material  is  said  to  exhibit  normality  if  the 

plastic  strain  increment  vector  dcP  is  parallel  to  the  normal  to  the  yield 
surface  regardless  of  the  location  of  the  interior  point  and  the  path  to  the 
surface.  If  a  system  comprised  of  this  material  and  loaded  by  an  arbitrary 
number  of  surface  tractions  is  observed  it  too  will  possess  a  yield  sur¬ 

face  in  a  force  space  having  a  dimension  equal  to  the  number  of  the  P^  .  De- 
noting  the  plastic  displacements  corresponding  to  the  P^  as  and  super¬ 

posing  the  axes  of  a  plastic  displacement  space  composed  of  the  uj^  on  the 
force  space  P^  ,  the  attainment  of  yielding  along  some  path  in  force  space 
will  again  be  accompanied  by  plastic  flow.  If  the  system  is  again  comprised 
of  a  material  exhibiting  normality,  the  plastic  displacement  increment  vector 
du^  will  also  be  parallel  to  the  normal  to  the  yield  surface  at  the  yield 

Y 

load  •.  Conversely  a  system  exhibiting  normality  in  force  space  is  com¬ 

prised  of  materials  exhibiting  normality  in  stress  space  provided  there  is 
no  frictional  dissipation  on  the  boundary  of  the  system. 

The  failure  envelope  for  sands  is  a  surface  definable  in  both  force  and 
stress  space.  A  lack  of  normality  to  this  surface  in  either  space  may  be 
explained  in  two  ways.  Either  the  failure  envelope  is  not  a  yield  surface 
or  the  material  itself  does  not  exhibit  normality.  Discussion  concerning 
the  first  of  these  concepts  is  given  in  Drucker  et  al  (9),  Jenike  and  Shield 
(18),  Drucker  (12)  and  Palmer  (21)  while  the  second  is  treated  in  Drucker  (8), 
Palmer  (21)  and  this  author  and  Paslay  (31).  The  purpose  of  this  investigation 
was  not  to  attempt  resolution  of  this  theoretical  question,  but  simply  to  study 
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the  character  of  the  flow  behavior  of  cohesionless  soils  at  failure. 

EXPERIMENTAL  PROCEDURE 

This  section  contains  a  discussion  of  the  experimental  procedure  followed 
in  this  investigation.  A  description  of  the  testing  device  used  and  of  any 
modifications  to  the  generally  accepted  procedures  for  this  device  is  included. 
Succeeding  sections  deal  with  the  testing  program,  the  results  obtained,  and 
the  analysis  of  the  data. 

Conceptually  the  direct  shear  device  provides  a  means  for  evaluating 
the  effect  of  shearing  stress  on  a  prescribed  plane  of  a  soil  as  a  function 
of  the  normal  stress  existing  on  that  plane  and  the  intrinsic  properties  of 
the  material.  In  the  late  thirties  and  probably  earlier,  enough  theoretical 
evidence  had  been  presented  to  indicate  that  this  simple  conceptual  picture 
of  what  the  results  of  direct  shear  testing  indicated,  possessed  only  an  accid¬ 
ental  connection  with  reality.  The  state  of  stress  and  strain  within  a  direct 
shear  test  specimen  is  exceedingly  complex  and  any  abstraction  of  the  normal 
and  tangential  surface  loads  and  their  associated  displacements  for  indicating 
pointwise  material  properties  can  not  be  substantiated.  Soil  researchers,  on 
acceptance  of  these  concepts,  turned  their  efforts  towards  development  of  an 
apparatus  theoretically  capable  of  predicting  point  material  properties  in  a 
finite  sized  specimen. 

In  this  study  only  the  forces  and  their  associated  displacements  acting 
on  the  direct  shear  specimen  are  recorded.  No  attempt  is  made  to  abstract 
these  results  to  the  stress-strain  level.  The  primary  object  of  the  study 
is  to  ascertain  the  relationship  between  the  direction  of  the  plastic  defor¬ 
mation  rate  vector  and  the  direction  of  the  normal  to  the  failure  envelope.  The 
inferences  as  to  pointwise  cohesionless  soil  properties  from  these  results  can 
only  be  qualitative.  Nevertheless  it  should  be  kept  in  mind  that  the  direct 
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shear  test  gives  a  serviceable  indication  of  the  magnitude  of  the  friction 
angle,  is  reflective  of  initial  void  ratio,  and,  perhaps,  most  importantly 
predicts  the  flow  in  force-displacement  terms  of  the  stress-strain  curves 
associated  with  triaxial  testing.  However,  only  further  testing  with  tri- 
axial  equipment  can  substantiate  the  inferences  made  from  the  results  present¬ 
ed  herein. 

Figs.  2,  3,  and  4  show  respectively  for  the  device  used  in  the  investi¬ 
gation;  a  vertical  cross  section  through  the  specimen  in  the  plane  of  the  ap¬ 
plied  tractions,  an  enlargement  of  this  cross  section,  and  a  plan  view  of  the 
apparatus.  The  various  parts  of  importance  in  the  ensuing  discussion  are  in¬ 
dicated  by  capital  letters  and  dimensions  required  for  calculation  by  letters 
in  lower  case.  The  views  shown  in  these  three  figures  represent  the  disposition 
of  the  system  prior  to  the  application  of  the  weight  (A)  and  the  horizontal 
thrust.  The  specific  device  used  in  the  study  was  a  Soiltest,  Inc.  Direct 
Shear  Test  Apparatus,  Model  D-110  equipped  with  a  Soiltest,  Inc.  Direct  Shear 
Box,  Model  D-183.  The  apparatus  has  been  modified  slightly  to  permit  recording 
of  the  necessary  data  required  in  this  investigation.  A  hand  operated  tangential 
thrust  mechanism  was  selected,  partially  for  better  control  over  incrementing 
the  tangential  load,  but  primarily  to  eliminate  the  effect  on  vertical  displace¬ 
ments  of  the  vibrations  of  a  previously  used  motorized  apparatus. 

Before  discussing  in  detail  the  particular  procedural  aspects  associated 
with  this  series  of  experiments,  a  brief  description  of  a  direct  shear  test 
is  included  to  provide  a  base  for  the  ensuing  discussion.  Reference  will  be 
made  to  the  lettered  parts  of  Figs.  2,  3,  and  4.  The  sample  (J)  is  placed 
at  some  initial  density  within  the  well  formed  by  the  upper  (H^)  and  lower  (H^) 
square  shear  box  rings  and  the  lower  loading  block  (D^)-  Parts  (H^)  and  (D2> 
are  rigidly  attached  to  the  shear  box  housing  (F)  and  this  combination  may  be 
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considered  fixed  with  respect  to  ground.  An  assembly,  which  during  the  test 
may  be  considered  rigid,  consisting  of  the  upper  loading  block  (D^),  the  yoke 
assembly  (C)  and  the  tilt  bar  (E)  is  then  applied.  In  the  normal  test  dial 
gages  (Ry)  and  (R^)  are  placed  in  the  positions  shown  to  record  respectively 
the  vertical  movement  of  the  tilt  bar  E  and  the  horizontal  displacement  of 
the  upper  square  shear  box  ring  (H^).  After  noting  the  initial  readings  of 
these  dials,  the  vertical  load  is  applied  to  the  weight  hanger  (B)  by  weights 
(A)  to  a  prescribed  value.  Subsequently,  after  pulling  pins  connecting  the 
two  shear  box  rings,  a  horizontal  thrust  is  applied  by  a  screw  mechanism  at 
the  right  of  Fig.  2.  The  magnitude  of  this  reactive  thrust  is  determined  using 
the  dial  gage  (Rp)  of  the  proving  ring  mounted  in  series  with  the  feed  screw. 
The  thrust  is  transmitted  from  the  proving  ring  to  the  reactor  shaft  and  thence 
to  the  shear  box  coupling  (G).  The  shear  box  coupling  is  rigidly  attached  to 
the  upper  square  box  ring  (F^).  The  reactive  thrust  drives  this  ring  to  the 
left  and  thus  deforms  the  specimen,  readings  of  the  gages  being  taken  at  inter¬ 
vals.  The  upper  loading  block  (D^)  remains  fixed  laterally  with  respect  to 
the  upper  shear  box  ring  (H^)  by  means  of  the  fins  protruding  from  the  loading 
block  into  the  sample .  During  the  loading  process ,  however ,  the  upper  loading 
block  (D^)  is  free  to  translate  vertically  or  to  rotate. 

The  sample  is  said  to  have  failed  whenever  the  tangential  thrust  required 
to  cause  a  increase  in  horizontal  deformation  begins  to  diminish.  This  conven¬ 
tional  definition  of  failure  in  direct  shear  fits  the  concept  of  system  inst¬ 
ability  found  in  Drucker  (11)  if  the  vertical  and  horizontal  loads  are  the 
only  tractions  doing  work  on  the  sample,  and  the  vertical  load  remains  of  con¬ 
stant  magnitude.  In  the  terminology  of  Drucker  (11)  a  system  is  said  to  be¬ 
come  unstable  whenever  the  increment  of  work  done  on  the  system  by  the  load 
increments  AP,.  acting  thru  their  corresponding  displacement  increments  Au. 
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becomes  negative.  Thus  if  readings  are  taken  at  intervals,  instability  of 

"t  h  *t  l*i 

the  system  occurs  in  the  interval  between  the  k  and  k  plus  one  readings 
if 

{P.  (k  +  1)  -  P.  (k) }  {u.  (k  +  1)  -  u. (k) }  <  0  (1) 

1  XI  1 

i.e. 

{Px  (k  +  1)  -  P1  (k) >  {u1  (k  +  1)  -  ux  (k)} 

+  {?2  (k  +  1)  -  ?2  (k) >  {u2  (k  +  1)  -  u2  (k)} 

+...+  {P  (k  +  1)  -  P  (k)}  {u  (k  +  1)  -  u  (k)}  <  0  (2) 

n  n  n  n 

for  a  system  subjected  to  n  surface  tractions  doing  work  on  the  specimen. 

For  the  specific  device  employed  in  this  investigation  the  stresses  in 
the  sample  tend  to  induce  a  rotation  of  the  upper  loading  block  (D2)  in  a 
clockwise  direction  about  an  axis  perpendicular  to  the  plane  of  the  diagram  of 
Fig.  2.  This  rotation  causes  an  appreciable  moment  on  the  upper  surface  of 
the  specimen  due  the  length  of  the  yoke  assembly  (C).  Evaluation  of  the  magni¬ 
tude  of  the  rotation  was  determined  by  inclusion  of  the  tilt  bar  (E)  and  the 
tilt  gage  (R^,)  through  a  calculation  to  be  mentioned  in  a  subsequent  paragraph. 
Although  this  was  the  only  additional  work  producing  force  included  in  this 
investigation  more  subtle  forms  exist.  For  the  present  study,  their  effect  on 
Inequality  (2)  and  on  the  lack  of  normality  was  felt  to  be  negligible.  Two  of 
these  tractions  which  might  be  of  consequence,  however,  are  the  possible  shear¬ 
ing  force  set  up  in  a  vertical  direction  between  the  left  hand  interior  wall 
of  the  upper  shear  box  ring  (H^)  and  the  specimen  due  to  a  relative  displace¬ 
ment  observed  near  failure  between  this  wall  and  the  specimen,  and  secondly  a 
moment  induced  by  rotation  of  the  upper  shear  box  ring  (H^)  about  an  axis  coin¬ 
cident  with  the  intersection  of  the  right  hand  exterior  wall  and  lower  surface 
of  the  ring  (H^).  As  the  determination  of  the  first  of  these  force  magnitude 
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is  somewhat  difficult  and  the  estimations  made  of  both  their  contributions 
being  of  small  magnitude,  the  assessment  of  these  two  effects  was  not 
included  in  the  testing  program. 

As  indicated  in  the  previous  paragraphs  three  force  variables  will  be 
considered  to  do  work  on  the  soil  system:  the  horizontal  or  tangential  traction 
T,  the  vertical  normal  load  N,  and  a  moment  M  induced  by  rotation  of  the 
upper  loading  block.  The  three  displacements  corresponding  to  these  forces 
are  respectively:  the  horizontal  displacement  6^  of  the  upper  shear  box  ring 
(H^),  the  vertical  displacement  of  the  midpoint  of  the  lower  surface  of 

the  upper  loading  block  (D.^),  and  the  rotation  \f>  of  this  surface  about  an 
axis  through  its  midpoint  and  normal  to  the  plane  of  the  tractions.  Fig.  5 
shows  a  unstressed  soil  specimen  and  the  positive  directions  of  these  forces 
and  their  associated  displacements.  In  Fig.  6  a  schematic  view  of  the  initial 
and  displaced  configurations  is  presented.  Using  Fig.  3,  5  and  6  as  references 
the  following  analytic  expressions  for  the  derived  variables  can  be  formulated. 
The  initial  tilt  of  the  sample  coincident  with  the  lower  edge  of  the  loading 
block,  designated  as  <|/q  is 

<|>0  =  sin"1  {aR  -  aL)/h}  (3) 

where  h  is  the  width  of  the  shear  box  ring  well  in  the  plane  of  Fig.  3  and 
the  dimensions  aR  and  a^  are  averages  of  two  depth  gage  readings  taken 
from  the  top  of  the  upper  yoke  assembly  bar  (C)  to  the  points  marked  by  the 
symbol  *  on  Fig.  4.  Using  these  same  depth  gage  readings  and  the  initial  tilt, 
the  initial  average  specimen  thickness,  tQ,  is  determined  from 

tQ  =  d1  -  [{d3  -  d2>  -{(aL  +  aR)/2}]/cos^Q)  (4) 

In  Eq.  (4)  the  dimensions  d^,  dg,  and  d^  are  defined  by  Fig.  3.  Knowledge 

of  the  specific  gravity  of  the  sand  G  ,  the  sample  weight  W  ,  the  planer 

s  s 


-11- 
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area  of  the  well  b  ,  the  initial  average  specimen  thickness,  tQ,  and 

the  current  vertical  displacement  5^  enables  a  computation  to  be  made 

of  the  initial  and  current  average  sample  void  ratio  e,  i.e. 

e  =  {G  v  h2  (t  -  6„)/W  }  -  1  (5) 

s  w  o  N  s 

The  vertical  normal  load  N  is  the  sum  of  the  weights  (A)  designated  by 
N  and  the  combined  weight  N  of  the  weight  hanger  (B),  the  yoke  assembly 

cl  O 

(C),  the  tilt  bar  (E),  and  the  upper  loading  block  (D^).  Thus 

N  =  N  +  N  (6) 

o  a 

Calculation  of  the  remaining  two  force  variables,  T  and  M,  and  the 

three  displacements  <$T,  6N  and  ^  requires  the  use  of  the  four  dial  gage 

readings  Rp,  Rp,  Ry  and  R^  .  Positive  readings  of  Rp  and  Rp  were  induced 

by  shortening  their  respective  stems  while  the  reverse  was  true  for  Ry  and 

R.J,  .  Denoting  the  difference  between  the  current  and  initial  readings  of 
these  four  gages  also  by  Fp ,  R^,  Ry  and  RT,  the  tangential  thrust  T  is 
determined  from 


T  =  k  *  Rd  (7) 

P  P 

kp  being  the  proving  ring  constant.  The  horizontal  displacement  <5^  as¬ 
sociated  with  the  thrust  T  is  simply  equal  to  Rp  .  The  current  tilt 
angle  i|»  (see  Fig.  6)  can  be  determined  from 


ij>c  =  tan-i{(aT  •  sin\|>o  +  RT  ~  RyVa^osi^}  (8) 

where  a^,  is  the  initial  distance  along  the  tilt  bar  (E)  between  the  center- 
lines  of  the  dial  gage  stems  Ry  and  R^  .  Hence  the  rotation  becomes 


ip  =  \p  -  ip  (9) 

c  o 

The  vertical  displacement  of  the  point  Q  of  Fig.  6,  denoted  as  6^,  and 
associated  with  the  vertical  normal  force  N  can  be  shown  to  be 


fi=  R..  +  6  *  ta#  +  d„{(cos^  -  cosi^  ) 

NvT  C3  C  O 

+  ta#  (si#  -  si#  )} 
o  c  o 


(10) 
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Determination  of  the  moment  M  caused  by  the  lateral  displacement  of  the 
weights  from  a  vertical  line  through  point  Q  is  given  by 

M  =  {N  (r  +  d  -  d  )  +  N  (r  +  d  -  d  )}  *  si n^  (11) 

o  c  z  6  a  z  o  c 

with  the  dimension  r  being  defined  by  Figure  3  and  (r  +  d  -  d  )  being 

G  Z.  o 

the  parallel  distance  from  point  Q  to  the  center  of  gravity  of  the  com¬ 
bination  of  parts  (B),  (C),  (E)  and  (D^). 

TESTING  PROGRAM 

The  particular  sand  selected  for  the  series  of  experiments  was  a  locally 
quarried  quartz  sand  commonly  used  in  the  making  of  concrete.  The  material 
was  screened  and  only  particles  passing  a  #20  mesh  sieve  (U.S.  Standard)  but 
retained  on  a  #25  mesh  sieve,  were  used  in  the  testing  program.  Under  the 
microscope,  the  grains  appeared  to  have  a  predominantly  angular  shape.  A 
determination  of  the  specific  gravity  of  the  sand  yielded  a  figure  of  2.65 
while  additional  tests  indicated  a  possible  void  ratio  range  of  0.662-0.965 
for  the  particles  making  up  the  samples. 

In  general  the  testing  program  involved  preparing  the  specimen  to  some 
prescribed  initial  configuration,  applying  the  normal  load  N  to  a  parti- 
cular  magnitude  and  finally  applying  the  tangential  thrust  T  in  increments, 
until  the  tangential  displacement  reached  a  value  between  0.25  -  0.38  centi¬ 
meters  (0.10  -  0.15  inches).  In  some  cases,  however,  this  procedure  was  al¬ 
tered  by  cycling  the  tangential  thrust,  changing  or  cycling  the  normal  load, 
or  a  combination  of  the  above.  Such  nonstandard  loading  was  included  parti¬ 
ally  to  ascertain  its  effect  on  failure  but  primarily  to  obtain  information 
regarding  the  inelastic  behavior  of  sands  within  the  failure  envelope.  This 
latter  subject  will  not  be  pursued  in  this  report. 

To  characterize  the  individual  tests,  a  group  symbol  has  been  assigned 


-13- 


to  each,  reflecting  the  initial  configuration  and  the  normal  load  N  at  fail¬ 
ure.  For  example,  5zd  indicates  that  the  normal  load  at  failure  is  of  the  5th 
magnitude  in  a  series  of  seven,  the  initial  tilt  z  is  in  the  zero  range,  and 
the  initial  average  void  ratio  d  is  in  the  dense  portion  of  the  spectrum. 
Amplifying  this  terminology  the  seven  final  normal  loads  are:  1  -  7.8  kg, 

2  -  11.8  kg,  3  -  19.8  kg,  4  -  35.8  kg,  5  -  51.8  kg,  6  -  67.8  kg  and  7  -  83.8  kg. 
The  initial  tilt  is  indicated  by  the  letter  symbols  p,  z,  and  n.  The 

letter  p  represents  +2.00°  z  ifi^z  +  0.68°,  while  z  designates  +0.67°*^* 

-  0.67°  and  n,  -  0.68°5  <1*  z  -  2.00°.  Lastly,  the  initial  average  void 

ratio  eQ  is  characterized  by  the  letters  d,  m,  and  1.  A  dense  initial  state, 
denoted  by  d,  requires  0.662  s  eQ  £  0.761,  a  medium  initial  state  m,  for 

0.762  $  e  £  0.865  and  a  loose  initial  state  1,  for  0.866  $  e  $  0.965.  Table  I 
o  o 

found  in  Appendix  III  depicts  this  grouping  of  the  107  test  numbers  making  up  the 
series  of  experiments. 

In  the  performance  of  an  actual  test,  readings  of  the  dial  gages  were 

taken  at  each  change  of  weight  (A)  and  at  each  increment  of  tangential  thrust. 

The  gages  were  read  in  a  systematic  manner  by  a  single  operator  and  thirty 

seconds  were  allowed  after  each  change  of  loading  to  permit  the  time  effect, 

attendant  to  granular  behavior,  to  subside.  The  size  of  each  tangential  thrust 

increment  was  in  the  neighborhood  of  3  kg  until  the  tangential  displacement 

increment  AT  reached  a  preset  value.  Control  of  the  thrust  mechanism  was 

then  shifted  to  the  displacement  gage,  R„  .  In  this  manner  a  well  spaced  set 

n 

of  values  was  obtained  for  plotting  purposes  and  the  work  increments,  defined 
by  Inequality  (2)  were  meaningful  for  locating  the  magnitude  of  the' variables 
associated  with  failure.  A  computer  was  utilized  in  reducing  the  raw  dial 
gage  readings  and  other  input  to  the  variables  defined  by  Inequality  (2)  and 
Eqs.  (3)  through  (11). 
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Plots  of  various  quantities  for  two  typical  tests  are  shown  as  Figs.  7, 
8,  9  and  10.  The  two  tests  selected  for  plotting  were  Nos.  452  and  472 
having  group  symbols  5zd  and  5zl  respectively.  Fig.  7  shows  projections  of 
the  three  dimensional  loading  histories  onto  the  T-N  and  T-M  force  planes 
while  Fig. 8  depicts  projections  the  corresponding  displacement  histories.  In 
Fig.  9  plots  of  the  tangential  thrust  versus  tangential  displacement  and  nor¬ 
mal  load  versus  normal  displacement  are  given.  Finally  Fig.  10  represents 
graphs  of  tangential  thrust  versus  the  configuration  variables  of  relative 
sample  density,  RD,  and  the  current  tilt  angle  \ .  The  relative  density 
is  defined  in  the  usual  manner  as 


RD  (%)  =  { (e  -  e)/  (e 


-  e  .  )} 
mm 


where  emax  and  e  ^  are  respectively  the  void  ratios  corresponding  the 
loosest  and  densest  states  (0.965  and  0.662). 

As  has  been  indicated  in  the  preceding  paragraphs  and  more  succinctly  by 
Table  I,  the  testing  program  reported  herein  involved  probing  the  soil  system 
under  investigation  over  a  wide  range  of  initial  conditions  and  loading  paths. 
A  large  number  of  tests  was  needed  in  order  to  provide  sufficient  data  for 
determining  a  mathematical  expression  which  would  adequately  represent  the 
failure  surface  in  a  hyperspace  of  the  three  force  variables  and  any  necess¬ 
ary  configuration  parameters.  The  results  of  the  experimental  program  are 
presented  in  the  following  section. 


TEST  RESULTS 

For  the  particular  soil  system  discussed  herein,  failure  of  the  specimen 

"til 

was  assumed  to  have  occurred  in  the  interval  between  the  k  and  the  kx  plus 
one  observations  whenever  Inequality  (2)  became  satisfied.  In  terms  of  the 
work  inducing  forces  acting  on  the  system  this  Inequality  may  be  rewritten  as 
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{T(k  +  1)  -  T(k)}*{6T(k  +  1)  -  <ST(k)}  +  {N  (k  +  1)  -  N(k)}*{6N(k  +  1)  -  6N(k) > 

+  {M(k  +  1)  -  M(k)}* {ij»(k  +  1)  -  i|»(k)}  $  0  (13) 

Once  this  inequality  is  satisfied  the  actual  forces  selected  as  those  acting 

"til 

at  failure  are  the  ones  associated  with  the  k  observation.  Since  the  magni¬ 
tudes  of  these  forces  change  only  slightly  near  failure  for  the  interval  sizes 
selected  the  maximum  percent  error  between  the  real  and  selected  force  values 
at  failure  was  estimated  to  be  2%  with  the  average  being  somewhat  less  than 
0.5%. 

The  values  determined  for  the  plastic  displacement  increments  associated 
with  failure  were  determined  in  the  following  manner  and  based  on  an  assumption 
regarding  the  elastic  displacements.  Previous  testing  had  indicated  that  on 
unloading  from  a  force  point  near  failure,  the  recoverable  or  elastic  displace¬ 
ments  were  zero  for  a  unloading  increment  corresponding  in  magnitude  to  the 

previous  loading  change.  Thus  the  displacements  occurring  on  loading  from  the 
"th  th 

k  1  to  the  k  plus  one  observation  were  totally  irrecoverable  or  plastic. 

Stated  in  another  way  this  assumption  implies  that  the  stiffness  moduli  as¬ 
sociated  with  reversible  behavior  are  three  to  four  orders  of  magnitude  larger 
than  the  moduli  associated  with  irreversible  behavior  at  failure  and  thus  for 
displacement  increments  of  0.01"  maximum,  the  elastic  displacements  would  re¬ 
main  undetected  in  dial  gage  readings  estimated  to  0.00001".  Based  on  this 
assumption  the  plastic  displacement  increments  were  simply  the  changes  A6^,, 
A6n,  Aip  occurring  in  the  interval  in  which  Inequality  (12)  became  satis¬ 
fied.  For  the  scale  of  Fig.  8  the  displacement  6^  and  the  rotation  i/»  ap¬ 
pear  to  be  fairly  smooth  functions  of  the  displacement  6^  near  failure  but 
on  a  larger  scale  they  are  somewhat  more  ragged.  Consequently  the  values  of 
A6n  and  actually  reported  represent  unweighted  smoothings  of  those 

portions  of  the  6^  -  6^  and  ip  —  6^,  curves  over  five  displacement  increments 
centered  on  the  increment  in  which  Inequality  (13)  was  initially  satisfied. 
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Table  II  of  Appendix  III  gives  the  configuration  data  of  all  the  tests 
reported  for  both  the  initial  and  failure  states.  The  test  number  indicates 
the  chronological  order  of  testing.  In  this  sequence,  the  missing  numbers 
are  indicative  of  tests  not  carried  to  failure.  The  sand  for  each  experi¬ 
ment  came  from  the  same  source,  was  not  reused,  and  all  tests  were  performed 
by  a  single  operator  over  a  three  month  period.  The  laboratory  temperature  is 
kept  at  21±  1°  C  with  the  relative  humidity  being  50±  10%.  All  samples  were 
tested  in  the  air  dried  condition.  The  forces  and  displacement  increments  at 
failure  are  listed  in  Table  III  of  Appendix  III.  In  the  subsequent  and  final 
section  of  this  paper  the  analysis  made  using  the  data  of  Tables  II  and  III  is 
presented  and  discussed. 


ANALYSIS  OF  TEST  DATA 

The  analysis  of  the  data  contained  in  Tables  II  and  III  is  essentially 
a  two  stage  process.  Initially  a  nest  of  failure  envelopes  in  the  force 
space  of  the  applied  tractions  must  be  defined  by  an  analytic  expression  from 
which  normals  may  be  determined.  Subsequently  a  systematic  comparison  between 
the  normals  and  the  corresponding  plastic  displacement  increment  vectors  must 
be  made.  In  accomplishing  this  analysis  process,  various  aspects  of  the 
statistical  method  are  employed.  These  statistical  techniques  facilitate  the 
handling  of  a  large  amount  of  data  and  more  importantly  are  useful  in  helping 
to  substantiate  inferences  made  from  the  results. 

In  the  first  stage  an  analytic  expression  was  sought  which  described  fail 
ure  of  the  system  in  terms  of  the  force  variables  and  the  configuration  para¬ 
meters.  The  technique  employed  to  determine  this  function  was  simply  to  fit  a 
surface  by  the  method  of  least  squares,  that  would  adequately  predict  the  tang 
ential  thrust  at  failure  in  terms  of  the  other  variables.  Pictorally  such  a 
function  would  yield  a  nest  of  failure  envelopes  fhe  force  space  of  the 
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applied  tractions.  A  single  surface  of  the  nest  would  be  isolated  by  holding 
constant  whatever  configuration  variables  were  necessary  in  defining  the 
failure  envelope  function.  As  the  sole  purpose  of  this  stage  was  to  deter¬ 
mine  an  analytic  function  from  which  normals  could  be  extracted  the  question 
as  to  whether  the  form  of  the  function  eventually  selected,  was  correct  theore¬ 
tically,  appeared  to  be  of  secondary  importance.  Certain  restrictions  of  a 
physical  nature  were,  however,  included  in  the  development  of  the  possible 
formulations.  The  selection  of  the  failure  tangential  thrust  as  the  dependent 
variable  is,  of  course,  somewhat  arbitrary  but  considering  the  nature  of  the 
system  tested,  the  choice  was  felt  to  be  logical. 

As  previously  mentioned,  the  independent  variables  were  taken  as  the 
failure  normal  force  and  moment,  and  additionally  those  system  configuration 
parameters  necessary  to  formulate  an  adequate  failure  function.  Such  config¬ 
uration  parameters  include  the  initial  density,  the  initial  tilt  angle,  and 
the  magnitudes  of  the  displacements  and  rotation  at  failure.  How  many  of  these 
parameters  and  what  combinations  of  them  would  be  required  to  formulate  an 
adequate  predictive  equation  for  the  tangential  failure  thrust  were  additional 
issues  to  be  resolved  in  this  first  stage  of  the  analysis. 

The  basic  form  of  the  expression  chosen  to  predict  the  tangential  thrust 
at  failure  was  an  infinite  power  series  in  the  independent  variables. 


n 

E 

i=l 


c. 

1 


(1) 


X. 

i 


n  n 
+  E  E 
i=l  j-1 


X.X. 
i  1 


+  . . . 


(14) 


In  Eq.  (14)  T^  is  the  predicted  value  of  tangential  thrust  at  failure,  X 
is  the  symbol  for  an  independent  variable  and  c  symbolizes  the  constants  to 
be  determined  by  the  method  of  least  squares.  The  actual  independent  variables 
required  for  an  adequate  predictive  equation  were  eventually  determined  to  be 
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the  vertical  normal  force  at  failure  N^,  the  moment  at  failure  M^,  and 
the  failure  void  ratio  e^.  .  In  reality  it  can  be  seen  from  Eq.  (11)  that 
M^.  is  a  function  of  and  the  current  tilt  angle  at  failure  .  Thus 

a  somewhat  more  basic  independent  variable  could  be  \p^,  but  since  the  re¬ 
quirements  of  the  second  stage  of  the  analysis  call  for  a  normal  vector  with  a 
component  along  the  M  axis,  M^.  has  been  taken  as  independent.  Some  of 
the  constants  c  can  be  set  to  zero  in  Eq.  (14)  if  it  is  assumed  that 
equals  zero  when  N^.  is  zero.  This  assumption  was  made  in  this  analysis 
and  thus  all  terms  in  the  expanded  version  of  Eq.  (14)  not  containing 
as  a  variable  are  absent.  The  lowest  powered  form  of  Eq.  (14)  appearing  to 
give  satisfactory  predictions  for  T^  was  the  third.  Rewriting  Eq.  (14)  in 
expanded  form  up  to  the  third  power  and  incorporating  the  aforementioned 
physical  assumption  concerning  N^_  yields 

T _  =  +  b„N  2  +  b-N_M_  +  b,  N-e.  +  bcN  3  +  b„N  _2M. 

fp  If  2  f  3  f  f  4  f  f  5  f  6ff 

*  Vf%  *  Wf2  +  VfVf +  bioNfef2  (15) 

where  the  b^  represent  the  constants  to  be  determined  from  the  method  of 
least  squares.  On  performing  the  numerical  analysis  the  actual  equation  ob¬ 
tained  from  the  input  of  Tables  II  and  III  is 

T _  =  -0.17749947  -0.01211425  N  2  +  0.00783852  N^M  +  5.20248222  N^e^ 

f p  f  f  f  f  f  f 

+0.00002449  Nf 3  -0.00001539  Nf2Mf  +0.00985996  Nf2ef  -0.00002242  NfMf2 
-0.00571806  NfMfef  -4.62136745  Nfef2  (16) 

The  standard  error  of  estimate  of  T^.  on  N^,  and  e^.,  symbolized 

by  s„  M  M  is  a  measure  of  the  scatter  about  the  regression  surface  of  T,. 

if*  fMfef 

on  N,.,  M,.  and  e,..  Surfaces  parallel  to  that  of  Eq.  (16)  ±s_  M  kg  from 
r  £  x 
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it  along  the  T  axis  should  contain  68%  of  all  the  sample  points  if  the 

randomness  is  normally  distributed, s„  „  M  is  defined  as 

r*  fref 


T£.N_M_e£  = 
f  f  f  f 


n 

Z 

i  =  1 


(Tfi  -  Vr 


(17) 


In  Eq.  (17)  n  is  the  number  of  tests,  the  other  symbols  being  previously 


defined.  For  the  data  of  this  series  s_  VT  was  found  to  be  1.8003  kg. 

TfWf 


A  somewhat  more  unbiased  estimate  of  the  standard  error  of  estimate  replaces 
the  denominator  of  Eq.  (16)  by  the  number  of  degrees  of  freedom.  As  ten  con¬ 
stants  exist  in  Eq.  (15)  and  nine  more  were  set  to  zero  by  the  physical  as¬ 
sumption,  the  number  of  degrees  of  freedom  is  (n-19).  Hence  the  unbiased  stand¬ 
ard  error  of  estimate  is  1.9851  kg. 

The  square  root  of  the  quotient  obtained  by  dividing  the  explained 
variation  of  by  the  total  variation  of  is  known  as  the  coefficient 

of  multiple  correlation,  R_  .In  this  expression  the  total  vari- 

f-  f  f®f 

ation  is  simply  the  variance  or  square  of  the  standard  deviation  of  T^.  while 

explained  variation  is  the  variance  of  T_  minus  the  square  of  s„  M  M 

2  f*  f  f  f 

Allowing  s„  to  stand  for  the  variance  of  the  coefficient  of  multiple 

f  t 

correlation  may  be  expressed  symbolically  as 


N  M  g 
f  f  f 


Wf 


(18) 


Performing  the  required  calculation  yielded  R  equal  to  0.9952. 

f.  f  f®f 

As  a  perfect  correlation  between  the  actual  failure  thrusts  and  those  pre¬ 
dicted  by  Eq.  (16)  would  be  unity,  it  was  felt  that  Eq.  (16)  adequately  re¬ 
presented  the  function  describing  the  nest  of  failure  envelopes. 
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The  dashed  lines  of  Figs.  11  through  16  show  intersections  of  failure 
surfaces  as  determined  by  Eq.  (16)  with  various  planes,  M  equal  a  constant 
or  N  equal  a  constant,  for  four  values  of  relative  density  (0%,  33%,  67% 
and  100%).  Also,  as  the  various  titles  of  the  figures  indicate,  the  infor¬ 
mation  contained  in  Tables  II  and  III  is  pictorally  displayed.  Cursory  ex¬ 
amination  of  the  figures  indicates  the  nature  of  the  agreement  between  Eq. 

(16)  and  the  actual  tangential  thrusts  at  failure.  Projections  of  the  plastic 
displacement  increment  vectors  onto  planes  where  N  is  constant  seem  to  be 
normal  to  the  failure  envelope  while  projections  of  these  same  vectors  onto 
planes  M  equal  a  constant  clearly  evidence  a  lack  of  normality. 

The  second  stage  of  the  analysis  was  directed  towards  systematically 
studying  the  relationship  between  the  direction  of  the  observed  plastic  de¬ 
formation  increment  vector  and  that  of  the  analytically  determined  normal  to 
the  failure  envelope.  For  a  given  test  the  directional  cosines  of  the  nor¬ 
mal  in  the  force  space  of  the  applied  tractions  were  obtained  by  differentiating 
Eq.  (16)  with  respect  to  N^.  and  .  The  partial  derivatives  3T^/3N^.  and 

3T^/3M^  thus  obtained  represent  slopes  of  the  projections  of  the  normal 
vector  onto  the  planes  M  =  0  and  N  =  0  respectively.  Inserting  the  parti¬ 
cular  values  of  N^,  M^,  and  e^  associated  with  the  test  in  question  into 
the  analytical  expressions  for  3T^/3N^  and  3T^/3M^.  yielded  numerical 
values  for  the  slopes  of  the  projections  from  which  the  directional  cosines  of 
the  normal  vector  could  be  calculated.  As  Eq.  (16)  represents  a  parametric 
nest  of  surfaces  in  force  space  best  fitting  the  data,  there  exists  the  some¬ 
what  obvious  observation  that  the  directional  cosines  calculated  by  the  pro¬ 


cedure  described  above  would  be  somewhat  different  from  those  obtained  by 


entering  Eq.  (16)  with  the  observed  values  T^,  and  .  This  latter 
method  would  yield  a  value  of  failure  void  ratio  different  from  the  observed 
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and  thus,  eventually,  modified  values  for  the  directional  cosines  of 
the  normal.  Actually  from  a  statistical  standpoint  a  more  elaborate  pro¬ 
cedure  would  be  necessary  to  obtain  a  less  biased  set  of  values  for  the 
directional  cosines  of  the  normals.  Such  a  procedure  would  involve  formul¬ 
ation  of  four  equations  of  the  type  exemplified  by  Eq.  (14)  with  each  express¬ 
ion  having  a  different  dependent  variable.  Relaxation  of  the  bias  introduced 
by  the  physical  assumption  which  set  certain  constants  to  zero  would  also  be 
necessary.  Manipulation  of  each  resulting  equation  would  collectively  yield 
four  sets  of  values  for  the  directional  cosines  of  any  particular  normal  that 
could  then  be  averaged.  Although  this  last  mentioned  method  is  theoretically 
more  sound,  the  close  correlation  of  Eq.  (16)  to  the  observed  data  as  express¬ 
ed  by  the  value  of  RT  strongly  suggests  that  the  simplified  techni- 

f-  f  f6f 

que  used  herein  to  determine  the  direction  of  the  normals  is  adequate. 

Having  the  directional  cosines  of  both  the  normal  and  flow  increment 
vectors  at  failure  for  each  test,  the  specific  purpose  of  this  investigation 
could  now  be  undertaken.  For  the  purpose  of  clarity  a  set  of  unit  vectors 
l,  ]  and  k  are  taken  to  act  respectively  along  the  M,  N,  and  T  axes  of  the 
force  space.  Denoting  the  directional  cosines  of  the  normal  to  the  failure 
envelope  to  the  M,  N  and  T  axes  by  l ,  m^  and  nn  respectively,  and  the 
directional  cosines  of  the  displacement  increment  vector  to  these  same  axes 
by  m^  and  n^  the  following  unit  vectors  may  be  defined  for  each  test 

result  as 

Plastic  Displacement  ^ 

Increment  Vector  =  =  l,t  +  m,^  +  n  Jc  (19. a) 

add 

Normal  Vector  =v  =  fcT  +  mT  +  n)c  (19.b) 

n  n  n 

Further  define  the  projections  of  these  vectors  onto  the  planes  M  equal  zero 
and  N  equal  zero  by  the  superscripts  NT  and  MT  respectively,  i.e. 
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A^1  =  m,  j  +  n,k  (20.  a) 

dJ  d 

A$MT  =  &  i  +  n,k  (20. b) 

d  d 

vNT  =  i  j  +  n  it  (20.c) 

nJ  n 

vMT  =  1  i  +  n  k  (20. d) 

n  n 

A  systematic  study  of  the  relationship  between  A?  and  v  for  each  test 


is  facilitated  by  calculation  of  the  following  angles. 

0A  =  cos'1  [(v*k)/(  | v  |  |k| )]  (21. a) 

6.,  =  cos"1  [(A$.k)/(|A5||k|)]  (21.  b) 

X5 

0C  =  cos  1  [(v*A^)/( |v | | A^| )]  (21. c) 

6d  =  cos"1  [(vMT*  k)/(  |vMT|  |k|)]  (21. d) 

0£  =  cos  1[(A^1  •  k)/(  | AS^"1 1  | k  1 )]  (21. e) 

6f  =  cos"1[(^T  •  A^T)/(|^T||t6MT|)]  (21‘f) 

e  =  cos"1  [(v1*1  .  k)/(|vNT||^!)3  (21. g) 

O 

e„  =  cos'1  C(A^NT  •  k)/( I a^nt 1 |k|)]  (21. h) 

n 

eT  =  cos'1  [(vNT  •  A^NT)/(|vNT|  |a^T|)]  (21. i) 


Fig.  17  illustrates  the  vectors  and  angles  defined  by  Eqs.  (19)  through  (21). 
Since  the  angles  given  by  Eqs.  (21.d)  through  (21.i)  are  restricted  to  lie 
in  the  planes  M  equal  zero  or  N  equal  zero  they  can  be  assigned  positive 
or  negative  values.  The  convention  adopted  is  based  on  the  following  equali¬ 
ties: 

0_  >  0  if  U  /n  )  >  0  (22. a) 

D  n  n 

6.,  >  0  if  (*,,/ n,)  >  0  (22. b) 

E  d  d 

0r  >  0  if  [U./n.)  “  (4/n  )3  >  0  (22. c) 

F  d  d  n  n 


0_  >  0  if  (m  /n  )  >  0 
G  n  n 


(22. d) 
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0  >  0  if  (m./n.)  >  0  (22. e) 

H  Q  Q 

0T  >  0  if  [(m,/n,)  -  (m  /n  )]  >  0  (22. f) 

i  da  n  n 

Ordering  of  angles  6.,  0O,  and  0  was  not  attempted  partially  because  such 

Ad  L 

a  process  is  somewhat  involved  but  primarily  because  sufficient  information 
was  obtained  by  simply  considering  their  absolute  magnitudes.  Calculation  of 
angles  0  and  0  is  somewhat  redundant  but  the  listing  of  these  values  en- 
ables  a  faster  visualization  of  the  behavior  of  the  normal  vector. 

Having  thus  defined  the  angular  relationships  expressed  by  Eqs.  (21), 
the  initial  part  of  this  second  phase  of  analysis  was  simply  to  compute  the 
numerical  values  of  these  quantities  for  each  of  the  tests  reported.  Once 
such  a  listing  is  determined  an  obvious  beginning  to  studying  the  character 
of  the  flow  is  simply  to  compute  the  means  and  standard  deviations  of  each 
of  the  angular  measures.  The  results  of  this  computation  are  given  in  Table 
IV  of  Appendix  III.  One  result  appears  to  stand  out  immediately.  The  mean 
of  the  angles  0p  denoted  as  0p,  is  very  close  to  zero.  Reference  to  Fig. 

'V* 

17  or  Eq.  (21.f)  discloses  that  such  a  value  for  0p  indicates  that  normality 
or  something  quite  close  to  it  exists  as  an  average  in  planes  where  N  equals 
a  constant.  A  second  result  of  somewhat  obvious  nature  is  that  normality  does 

<V 

not  exist  in  the  planes  of  constant  M  .  This  is  shown  by  the  value  of  0^, 
which  has  a  magnitude  of  -26.72°.  The  value  of  of  26.76°  reflects  both 

the  magnitudes  of  0^  and  0p  .  Geometrically  0^  represents  the  average 
co-latitude  of  the  flow  vector  from  a  normal  vector  coincident  with  a  polar 
axis.  The  importance  of  the  near  normality  result  indicated  by  0  is  es¬ 
sentially  to  place  all  of  the  variance  from  normality  indicated  by  the  magni- 
tude  of  into  planes  of  constant  M  .  For  the  system  studied  in  this  in¬ 

vestigation,  then,  lack  of  normality  can  be  studied  by  simply  concentrating 
on  the  projections  of  the  normal  and  flow  vectors  onto  the  plane  of  zero 
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moment.  Before  this  can  be  done,  however,  one  final  check  should  be  made  on 

t\f 

the  distribution  of  Qp  •  This  check  involves  determining  if  the  various 

*\i 

values  of  0p  about  0^,  vary  randomly,  or  whether  functional  relationships 
exist  between  the  0p  and  the  various  force  or  system  parameters  existing 
at  failure.  Two  parameters,  and  RD^,  were  selected  for  this  check. 

The  results  of  these  plots  are  shown  as  Figs.  18  and  19.  The  first  of  these 
figures  depicts  a  plot  of  the  various  0p  versus  their  corresponding  RD^. 
and  indicates  no  functional  relationship  existing.  The  second  figure,  however, 
does  appear  to  show  some  slight  dependence  of  the  0p  on  i .  This  depend¬ 
ence  indicates  that  the  divergence  from  normality  in  planes  of  constant  N 
may  increase  with  increasing  positiveness  of  at  failure.  Accompanying 

this  trend  and  making  it  somewhat  more  difficult  to  interpret  is  a  decided 
increase  in  the  deviations  about  the  mean  at  both  ends  of  the  investigated 
ipf  spectrum.  As  Fig.  19  does  not  appear  to  conclusively  demonstrate  depend¬ 
ence  of  0p  on  ipf  it  will  be  assumed  that  near  normality  with  random  vari¬ 
ance  exists  in  planes  of  constant  normal  force. 

The  character  of  the  definite  lack  of  normality  in  planes  of  constant 

•  •  «  'V 

moment  as  indicated  by  the  magnitude  of  0^.  may  also  be  investigated  by  the 
techniques  of  the  previous  paragraph.  As  previous  effort  in  this  area  both 
in  triaxial  and  direct  shear  testing  has  given  rise  to  the  phrase  "constant 
volume  flow,"  a  phrase  which  implies  that  the  rate  of  volume  change  at  fail¬ 
ure  is  zero  for  certain  relative  densities,  it  is  of  interest  to  also  study 
the  behavior  of  the  angle  between  the  projection  of  the  displacement  increment 
vector  onto  the  M  equal  zero  plane  and  the  T  axis,  i.e.  0  .  The  ob- 

vious  parameter  of  importance  here  is  the  final  relative  density  and  plots 
of  0^  and  0^  versus  this  parameter  are  shown  in  Fig.  20.  As  is  evident 
from  the  graph  both  0  and  0„  at  failure  are  dependent  on  the  relative 

1  n 
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density.  Assuming  a  linear  trend  for  both  plots  the  following  regression 
lines  were  determined. 


0U  =  2.28  -  0.2582  RD£  (23. a) 

n  r 

0X  =  -33.36  +  0.1184  RDf  (23. b) 

These  had  correlation  coefficients  as  determined  by  an  appropriately  modified 
form  of  Eq.  (18)  of  0.905  and  0.792  respectively.  The  ninety-five  %  confid¬ 
ence  limits  on  the  slope  of  Eqs.  (23. a)  and  (23. b)  were  ±0.0817  and  ±0.1179. 

A  closer  look  at  Figure  20  indicates  that  a  parabolic  curve  might  be  a  better 
trend  for  0„  versus  RD-  but  that  little  would  be  gained  by  use  of  such  a 

il  i 

form  for  0^  versus  RD^  .  Calculations  were  not,  however,  carried  out  to 
substantiate  these  observations. 

DISCUSSION  AND  SUMMARY 

For  the  system  investigated  normality  to  the  failure  surface  appears 
to  exist  for  projections  of  the  flow  increment  vector  onto  planes  of  constant 
normal  load.  This  is  indicated  by  the  nearly  zero  mean  value  of  0^,  and  the 
absence  of  any  conclusive  functional  dependence  on  the  various  force  and 
configuration  failure  parameters.  On  projecting  the  flow  increment  vector 
onto  planes  of  constant  moment  a  decided  lack  of  normality  is  seen.  As  a 
first  approximation  the  magnitude  of  the  angular  deviation  may  be  regarded 
as  a  constant  with  a  mean  value  of  -26.72°  and  a  standard  deviation  of  ±3.21°. 
Closer  examination  reveals,  however,  that  a  functional  relationship  exists 
between  0^  and  RD^.  .  This  can  be  seen  in  Fig.  20  or  determined  from 
Eq.  (23.b).  The  slope  term  in  Eq.  (23,b)  indicates  that  the  probable 
angular  deviation  from  normality  has  a  range  of  values  confined  to  a  fan 
having  a  central  angle  of  11.32°.  Most  natural  sands  have  relative  densities 
within  the  extremes  of  0  and  100%  and  thus  the  plausibility  of  using  a  con- 
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stant  value  (say  -  27.44°  at  RDf  =  50%)  to  describe  the  deviation 
from  normality  is  not  unrealistic. 

Before  attempting  to  draw  any  inferences  from  these  results  regarding 
the  possible  pointwise  behavior  of  a  cohesionless  soil,  two  important  limi¬ 
tations  of  the  testing  program  described  in  this  paper  should  be  discussed. 

The  first  of  these  is  that  the  final  loading  path  to  failure  in  all  tests 
was  confined  to  planes  of  constant  normal  load.  Consequently,  whether  or  not 
the  various  results  are  path  independent  has  not  been  shown.  The  second 
limitation  is  simply  that  since  only  one  specific  sand  was  tested  in  a 
single  device  there  can  be  no  assurance  that  the  results  can  be  applied 
to  a  wide  spectrum  of  granular  media.  Beside  these  limitations  it  perhaps 
should  be  repeated  that  inferences  drawn  from  observations  on  a  system  can 
only  be  qualitative  in  their  applicability  towards  the  formulation  of  stress- 
strain  relations. 

The  development  of  these  inferences  may  be  made  clearer  by  reference  to 
Fig.  21.  The  principal  stress  space  ,  o 2,  03  shown  in  Fig.  21. a  assumes 
compressive  stresses  to  be  positive.  Hydrostatic  pressures  are  located  on 
the  space  diagonal  where  =  a2  =  CT 3  *  On  superposing  the  axes  of  a 
principal  plastic  strain  space  e^,  £2^,  £3^  on  their  counterparts  in 
principal  stress  space,  plastic  volume  increases  have  the  same  direction  as 
that  of  the  negative  space  diagonal.  The  it-plane  in  Fig.  21. a  is  perpendi¬ 
cular  to  the  space  diagonal.  A  stress  or  plastic  strain  increment  vector  lying 
in  the  tt  plane  has  no  hydrostatic  or  dilatation  rate  component.  The  simplest 
analytical  representations  of  the  failure  surface  in  principal  stress  space 
for  cohesionless  soils  are  cones  or  pyramids  with  hexagonal  bases.  The  apexes 
of  these  surfaces  are  located  at  the  origin  of  coordinates  and  they  open  out¬ 
ward  with  increasing  hydrostatic  pressure  maintaining  a  minimum  of  three-fold 


-27- 


symmetry  about  the  space  diagonal.  The  rate  at  which  these  surfaces  diverge 
from  the  space  diagonal  is  a  function  of  the  current  relative  density  of  the 
material.  Intersections  of  three  of  these  surfaces  with  the  ir -plane  for  a 
specific  relative  density  are  shown  in  Fig.  21. b. 

The  implications  of  the  results  of  this  become  clearer  if  the  normal 
force  axis  (N)  of  the  force  space  of  applied  tractions  acting  on  the  direct 
shear  specimen  is  superposed  on  the  hydrostatic  axis  of  the  principal  stress 
space.  That  this  can  be  done  is  based  on  the  fact  that  the  system  volume 
change  is  reflected  only  in  the  displacement  6^  .  The  first  implication 

is  that  since  normality  is  preserved  in  planes  of  constant  normal  force  it 
is  preserved  in  the  ir-plane  of  principal  stress  space.  The  second  impli¬ 
cation  is  that  the  lack  of  normality  in  planes  perpendicular  to  the  ir-plane 
may  be  of  approximately  constant  angular  magnitude. 
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=  weights  applied  to  produce  normal  load  (see  Fig.  2); 

=  average  of  left  side  depth  gage  readings  (see  Figs.  3,  4,  and  6); 

=  average  of  right  side  depth  gage  readings  (see  Figs.  3,  4,  and  6); 

=  distance  along  tilt  bar  from  the  intersections  of  the  centerlines 
of  the  vertical  and  tilt  dial  gages  with  top  surface  of  tilt  bar 
(see  Fig.  3); 

=  weight  hanger,  direct  shear  device  (see  Fig.  2); 

=  symbol  for  constants  to  be  determined  by  method  of  least  squares 
(i  =  i,  2,  3  ...); 

=  yoke  assembly,  direct  shear  device  (see  Figs.  2,  3,  and  5); 

=  general  symbol  for  constants  in  Taylor  series; 

=  upper  loading  block,  direct  shear  device  (see  Figs.  2,  3,  and  5); 

=  lower  loading  block,  direct  shear  device  (see  Figs.  2  and  3); 

=  grouping  symbol  for  specimens  having  an  initial  relative  density 
between  66.7%  and  100.0%. 

=  generalized  plastic  displacement  increment  vector; 

=  plastic  strain  increment  vector  (see  Fig.  1); 

=  depth  dimension  from  top  surface  of  upper  shear  box  ring  (H^)  to 
upper  surface  of  lower  loading  block  (D2)  (see  Fig.  3); 

=  measured  thickness  of  tilt  bar  (E)  (see  Fig.  3); 

=  combined  measured  thickness  of  top  bar  of  yoke  assembly  (C), 
upper  loading  block  ( )  and  tilt  bar  (E)  (see  Fig.  3); 

=  tilt  bar,  direct  shear  device  (see  Figs.  2,  3,  and  4); 

=  average  void  ratio  of  specimen; 

=  initial  average  void  ratio  of  specimen; 

=  void  ratio  of  specimen  at  failure; 

=  loosest  void  ratio  obtainable  for  sand  used  in  experiments 
=  0.965; 

=  densest  void  ratio  obtainable  for  sand  used  in  experiments 
=  0.662; 

=  shear  box  housing,  direct  shear  device  (see  Figs.  2,  3,  and  4); 

=  function  of; 

=  shear  box  coupling,  direct  shear  device  (see  Figs.  2  and  3); 

=  specific  gravity  of  sand  used  in  experiments  =  2.65; 

=  upper  shear  box  ring,  direct  shear  device  (see  Figs.  2,  3,  and  4); 

=  lower  shear  box  ring,  direct  shear  device  (see  Figs.  2  and  3); 

=  interior  horizontal  dimension  of  square  shear  box  ring  (see 
Figs.  2  and  5); 

=  subscript  (i  =  1,  2,  3  ...); 

=  unit  vector  along  moment  axis  of  general  force  space; 

=  specimen,  direct  shear  device  (see  Figs.  2  and  3); 
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APPENDIX  II.  -  NOTATION  (continued) 
subscript  (j  =  1,  2,  3  . ..); 

unit  vector  along  normal  load  axis  of  general  force  space; 
right  hand  side  of  yield  function  equation  (see  Fig.  1); 
index  subscript  denoting  k**1  observation  of  dial  gages  during 
test; 

unit  vector  along  tangential  force  axis  of  general  force  space; 
proving  ring  constant; 

grouping  symbol  for  specimens  having  an  initial  relative  density 
between  0.0%  and  33.3%; 

direction  cosine  between  displacement  increment  vector  at 
failure  and  moment  axis; 

direction  cosine  between  the  normal  vector  to  the  failure  envelope 
and  the  moment  axis; 

moment  acting  on  specimen  induced  by  rotation  of  specimen  upper 

surface  (see  Fig.  5); 

moment  acting  on  specimen  at  failure; 

grouping  symbol  for  specimens  having  an  initial  relative  density 
between  33.3%  and  66.7%; 

direction  cosine  between  displacement  increment  vector  at  failure 
and  normal  load  axis; 

direction  cosine  between  the  normal  vector  to  the  failure  envelope 
and  the  normal  load  axis; 

normal  or  vertical  load  acting  on  specimen  (see  Fig.  5); 
that  portion  of  the  normal  load  on  the  specimen  due  to  the 
weights  (A); 

normal  load  acting  on  specimen  at  failure; 

that  portion  of  the  normal  load  on  the  specimen;  due  to  the 
combined  weights  of  the  direct  shear  device  parts  (B),  (C), 

(D1)  and  (E); 

grouping  symbol  for  specimens  having  an  initial  upper  surface 
tilt  angle  between  -0.68°  and  -2.00°; 

direction  cosine  between  displacement  increment  vector  at  failure 
and  the  tangential  force  axis; 

direction  cosine  between  the  normal  vector  to  the  failure  envelope 
and  the  tangential  force  axis; 
generalized  force  variable; 

yield  load  in  force  space; 

grouping  symbol  for  specimens  having  an  initial  upper  surface 
tilt  angle  between  +0.68°  and  +2.00°; 

symbol  for  tangential  displacement  dial  gage  and  the  reading 
of  this  gage  (see  Figs.  2  and  3); 

symbol  for  proving  ring  dial  gage  and  the  reading  of  this 
gage  (see  Figs.  2  and  4); 

symbol  for  tilt  dial  gage  and  the  reading  of  this  gage  (see 
Figs.  2  and  3); 

coefficient  of  multiple  correlation  of  the  tangential  force 
at  failure  on  the  normal  force,  moment  and  void  ratio  at  failure; 
symbol  for  vertical  displacement  dial  gage  and  the  reading  of 
this  gage  (see  Figs.  2  and  3); 
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=  average  relative  density  of  specimen; 

=  average  relative  density  of  specimen  at  failure; 

=  initial  average  relative  density  of  specimen; 

=  distance  along  centerlines  of  yoke  assembly  bars  (C)  from  lower 
surface  of  tilt  bar  (E)  to  contact  point  of  parts  (B)  and  (C); 

=  distance  along  centerlines  of  yoke  assembly  bars  (C)  from  lower 
surface  of  tilt  bar  (E)  to  center  of  gravity  of  parts  (B),  (C), 

(D)  and  (E); 

=  standard  deviation  of  the  variable  T_; 

f 

=  standard  error  of  estimate  of  the  tangential  force  at  failure 
on  the  normal  load,  moment  and  void  ratio  at  failure; 

=  tangential  or  horizontal  force  acting  on  specimen  (see  Fig.  5); 

=  tangential  force  on  specimen  at  failure; 

=  tangential  force  at  failure  predicted  by  regression  equation; 

=  generalized  displacement; 

=  generalized  plastic  displacement; 

=  weight  of  specimen; 

=  independent  variable  in  a  regression  analysis; 

=  dependent  variable  in  a  regression  analysis; 

=  grouping  symbol  for  specimens  having  an  initial  upper  surface 
tilt  angle  between  +0.67°  and  -0.67°; 

=  unit  weight  of  water; 

=  displacement  increment  in  direction  of  normal  force  at  failure ; 

=  displacement  increment  in  direction  of  tangential  force  at  failure; 

=  unit  plastic  displacement  increment  vector  at  failure  (see  Fig.  17); 

=  projection  of  unit  plastic  displacement  increment  vector  at 
failure  on  plane  N  =  0; 

=  projection  of  unit  plastic  displacement  increment  vector  at 
failure  on  plane  M  =  0; 

=  rotation  increment  in  direction  of  moment  at  failure; 

=  displacement  in  direction  of  normal  force  N  (see  Fig.  5); 

=  displacement  in  direction  of  tangential  force  T  (see  Fig.  5); 
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principal  plastic  strains; 

angle  between  normal  vector  v  and  T  axis  (see 


angle  between  displacement  increment  vector  a£ 
(see  Fig.  17); 

angle  between  v  and  A<5  (see  Fig.  17); 

“^MT 

angle  between  v  and  T  axis  (see  Fig.  17); 


angle  between 
angle  between 


a*mt  and  T  axis  (see  Fig.  17); 
vMT  and  A*MT  (see  Fig.  17); 
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and  T  axis  (see  Fig.  17); 
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APPENDIX  II.  -  NOTATION  (concluded) 


angle  between  a*nt  and  T  axis  (see  Fig.  17); 
angle  between  and  a*nt  (see  Fig.  17); 

unit  normal  vector  to  the  failure  envelope  (see  Fig.  17); 

projection  of  v  onto  plane  N  =  0; 

projection  of  v  onto  plane  M  =  0; 
general  state  of  stress; 

general  state  of  stress  causing  yield; 
principal  stresses; 

angle  of  internal  friction  at  failure  as  determined  in  direct 
shear  test; 

angle  of  internal  friction  at  failure  in  a  direct  shear  test 

showing  flow  at  constant  volume; 

rotation  in  direction  of  moment  M  (see  Fig.  5); 

current  angular  deflection  of  specimen  upper  surface  from 

horizontal  or  current  angle  of  tilt  (see  Fig.  6); 

initial  angular  deflection  of  specimen  upper  surface  from 

horizontal  or  initial  angle  of  tilt  (see  Figs.  5  and  6); 

angle  of  tilt  at  failure. 


APPENDIX  III.  -  TABLE  I 


GROUPING  OF  TESTS  BASED  ON  NORMAL  LOAD  AT 
FAILURE  AND  INITIAL  CONFIGURATION 


*  tests  which  included  internal  cycling  or  a  non-typical  loading  path 


APPENDIX  III 


TABLE  II 


INITIAL  AND  FINAL  CONFIGURATION  DATA 


Test  Group  Initial  Conditions  Conditions  at  Failure 

Num-  Sym-  ~  ~~~ 


ber 

bol 

Void 

Relative 

Ratio  e 

o 

Density 

RD  (%) 

0 

401 

3zd 

0.728 

78 

402 

4zd 

0.722 

80 

403 

5zd 

0.743 

73 

404 

6zm 

0.769 

65 

405 

7zd 

0.750 

71 

421* 

4zm 

0.766 

66 

422* 

4zd 

0.729 

78 

423* 

4zd 

0.730 

77 

424* 

4zd 

0.756 

69 

425* 

6zd 

0.728 

78 

426* 

6zd 

0.753 

70 

427* 

5zd 

0.750 

71 

428* 

5zd 

0.750 

71 

429* 

7zd 

0.714 

83 

430* 

7zd 

0.715 

82 

431 

4zl 

0.877 

29 

432 

6zm 

0.854 

37 

433* 

4zl 

0.923 

14 

434* 

6zl 

0.924 

13 

435* 

6zd 

0.741 

74 

436* 

6zd 

0.729 

78 

437* 

6zd 

0.734 

76 

438* 

6zd 

0.727 

78 

439* 

6zd 

0.719 

81 

440* 

5zd 

0.729 

78 

441* 

7zd 

0.723 

80 

442 

3pd 

0.693 

90 

443 

4pd 

0.723 

80 

444 

5pd 

0.675 

96 

445 

6pd 

0.687 

92 

446 

6nd 

0.690 

91 

447 

5nd 

0.685 

92 

448 

4nd 

0.681 

94 

449 

3nd 

0.669 

98 

450 

3zd 

0.726 

79 

451 

4zd 

0.742 

73 

Tilt 

Void 

Relative 

Tilt 

^ngle 
(deg.  ) 

Ratio  e^. 

Density 
RDf  (%) 

Angle  ip 
(deg.  ) 

0.04 

0.752 

70 

-0.38 

0.23 

0.740 

74 

-0.37 

1 

o 

• 

0.750 

71 

-0.80 

-0.30 

0.774 

63 

-1.19 

-0.11 

0.762 

67 

-0.92 

0.05 

0.783 

60 

-0.94 

-0.16 

0.759 

68 

CO 

CO 

• 

o 

0.07 

0.755 

69 

-0.59 

-0.30 

0.782 

60 

-1.07 

0.01 

0.748 

71 

-0.58 

-0.34 

0.764 

66 

-1.39 

-0.33 

0.772 

64 

-0.89 

-0.09 

0.770 

64 

-0.70 

-0.44 

0.735 

76 

-0.94 

0.16 

0.733 

76 

-0.29 

-0.05 

0.887 

26 

-0.49 

-0.10 

0.855 

36 

-0.51 

-0.07 

0.908 

19 

-0.72 

CO 

CM 

• 

o 

1 

0.903 

20 

-0.61 

-0.48 

0.764 

66 

-1.00 

-0.44 

0.754 

70 

-0.92 

-0. 34 

0.750 

71 

-0.83 

-0.21 

0.747 

72 

-0.79 

0.25 

0.742 

73 

-0.30 

0.09 

0.753 

70 

-0.63 

-0.25 

0.747 

72 

-0.81 

1.30 

0.726 

79 

0.71 

1.53 

0.744 

73 

0.55 

1.74 

0.691 

90 

1.07 

1.62 

0.698 

88 

1.01 

-1.59 

0.706 

85 

-1.72 

-1.70 

0.704 

86 

-1.80 

-1.64 

0.699 

88 

-1.70 

-1.63 

0.693 

90 

-1.72 

-0.24 

0.743 

73 

-0.64 

0.11 

0.754 

70 

-0.66 
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*  tests  which  included  internal  cycling  or  non-typical  loading  path 
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INITIAL  AND  FINAL  CONFIGURATION  DATA 


Test  Group  Initial  Conditions  Conditions  at  Failure 


Num¬ 

Sym¬ 

ber 

bol 

Void 

Relative 

Tilt 

Void 

Relative 

Tilt 

Ratio  e 

o 

Density 

RD  (%) 
o 

Angle  'P 
(deg.  )° 

Ratio  e^ 

Density 
RDf  (%) 

Angle  'J' 
(deg.  ) 

452 

5zd 

0.712 

83 

-0.16 

0.734 

76 

-0.52 

453 

6zd 

0.716 

82 

-0.16 

0.721 

80 

-0.57 

454 

3nd 

0.663 

100 

-1.67 

0.684 

93 

-1.77 

455 

4nd 

0.688 

91 

-1.57 

0.706 

85 

-1.69 

456 

5nd 

0.703 

86 

-1.47 

0.720 

81 

-1.62 

457 

6nd 

0.681 

94 

-1.78 

0.687 

92 

-1.77 

458 

3pd 

0.702 

87 

1.37 

0.720 

81 

0.77 

459 

4pd 

0.678 

95 

1.33 

0.697 

88 

0.90 

460 

5pd 

0.716 

82 

1.61 

0.727 

78 

0.90 

461 

6pd 

0.722 

80 

1.52 

0.736 

75 

0.64 

462 

lzd 

0.729 

78 

-0.30 

0.746 

72 

-0.64 

463 

2zd 

0.718 

81 

-0.04 

0.733 

76 

-0.44 

464 

3pm 

0.850 

38 

1.16 

0.858 

35 

0.40 

465 

4pm 

0.850 

38 

1.45 

0.854 

37 

0.57 

466 

5pl 

0.865 

33 

1.69 

0.858 

35 

1.14 

467 

6pl 

0.965 

0 

1.77 

0.916 

16 

1.33 

468 

3zl 

0.895 

23 

-0.12 

0.885 

26 

-0.51 

469 

4zl 

0.903 

20 

-0.29 

0.894 

23 

-0.60 

470 

3zl 

0.874 

30 

-0.02 

0.884 

27 

-0.27 

471 

4zl 

0.902 

21 

0.04 

0.898 

22 

-0.21 

472 

5zl 

0.900 

21 

0.07 

0.894 

23 

-0.34 

473 

6zl 

0.922 

14 

-0.15 

0.896 

23 

-0.56 

474 

3pm 

0.807 

52 

1.73 

0.823 

47 

1.30 

475 

4pm 

0.809 

51 

1.73 

0.826 

46 

1.04 

476 

5pl 

0.929 

12 

1.67 

0.890 

25 

0.83 

477 

6pm 

0.797 

55 

1.73 

0.808 

52 

0.72 

478 

3nm 

0.858 

35 

-1.54 

0.879 

28 

-1.50 

479 

4nd 

0.755 

69 

-1.79 

0.775 

63 

-1.86 

480 

5nl 

0.866 

33 

-1.78 

0.866 

33 

-1.73 

481 

6nm 

0.838 

42 

-1.80 

0.843 

40 

1.71 

482 

2zm 

0.815 

49 

-0.04 

0.828 

45 

-0.55 

483 

3zm 

0.841 

41 

-0.01 

0.840 

41 

-0.50 

484 

4zm 

0.826 

46 

-0.12 

0.827 

45 

-0.70 

485 

5zm 

0.840 

41 

-0.25 

0.834 

43 

-0.87 

486 

6zm 

0.796 

56 

-0.10 

0.806 

52 

-0.77 

490 

6nd 

0.728 

78 

-1.43 

0.751 

70 

-1.43 

APPENDIX  III.  -  TABLE  II  (concluded) 

INITIAL  AND  FINAL  CONFIGURATION  DATA 

Test  Group  Initial  Conditions  Conditions  at  Failure 


Num¬ 

Sym¬ 

ber 

bol 

Void 

Relative 

Tilt 

Void 

Relative 

Tilt 

Ratio  e 

Q 

Density 

Angle 

Ratio  e 

Density 

Angle  t 

RD  (%) 
o 

(deg.  ) 

X 

RDf  (%) 

(deg. ) 

491 

6zd 

0.719 

81 

-0.26 

0.742 

73 

-0.71 

492 

6pd 

0.750 

71 

1.06 

0.762 

67 

-0.52 

493 

5nd 

0.741 

74 

-0.94 

0.762 

67 

-1.27 

494 

5zm 

0.774 

63 

-0.19 

0.779 

61 

-0.99 

495 

5pd 

0.739 

74 

1.20 

0.750 

71 

0.17 

496 

4zd 

0.729 

78 

-0.16 

0.752 

70 

-0.57 

497 

4zm 

0.777 

62 

-0.09 

0.790 

58 

-1.05 

498 

4pd 

0.754 

70 

1.39 

0.767 

65 

0.61 

499 

3nm 

0.774 

63 

-1.07 

0.786 

59 

-1.40 

500 

3zd 

0.745 

72 

0.11 

0.762 

67 

-0.34 

501 

3pd 

0.753 

70 

0.97 

0.776 

62 

-0.04 

502 

6nm 

0.849 

38 

-1.27 

0.846 

39 

-1.26 

503 

6zm 

0.859 

35 

0.20 

0.859 

35 

-0.15 

504 

6pl 

0.868 

32 

0.95 

0.861 

34 

0.32 

505 

5nl 

0.885 

26 

-1.40 

0.880 

28 

-1.42 

506 

5zl 

0.900 

21 

-0.20 

0.888 

25 

-0.78 

507 

5pl 

0.887 

26 

1.33 

0.876 

29 

0.65 

508 

4nl 

0.876 

29 

-1.42 

0.873 

30 

-1.35 

509 

4zl 

0.893 

24 

-0.22 

0.888 

25 

-0.63 

510 

4pl 

0.919 

15 

1.15 

0.900 

21 

0.26 

511 

3nl 

0.876 

29 

-1.35 

0.878 

29 

-1.33 

512 

3zl 

0.877 

29 

-0.12 

0.883 

27 

-0.50 

513 

3pl 

0.890 

25 

0.89 

0.889 

25 

0.27 

514 

6zm 

0.825 

46 

0.11 

0.818 

48 

-0.42 

515 

5zm 

0.836 

42 

0.07 

0.833 

43 

-0.70 

516 

4zm 

0.812 

50 

0.13 

0.815 

49 

-0.46 

517 

3zm 

0.808 

52 

0.17 

0.823 

47 

-0.52 

518 

6nm 

0.821 

47 

-1.23 

0.  819 

48 

-0.89 

519 

5nm 

0.814 

50 

-1.11 

0.818 

48 

-1.57 

520 

4nm 

0.823 

47 

-1.29 

0.834 

43 

-1.42 

521 

3nm 

0.826 

46 

-1.08 

0.836 

42 

-1.13 

522 

6pm 

0.833 

43 

1.17 

0.821 

47 

-0.36 

523 

5pm 

0.820 

48 

0.85 

0.823 

47 

-0.65 

524 

4pm 

0.808 

52 

0.94 

0.809 

51 

0.08 

525 

3pm 

0.796 

56 

1.67 

0.810 

51 

0.19 

APPENDIX  III.  -  TABLE  III 


DISPLACEMENT  INCREMENTS  AND  FORCES  AT  FAILURE 


Test 

Group  Normal 

Shear 

Moment 

Vertical 

Horizontal 

Num¬ 

Sym¬ 

Force 

Force 

Displacement 

Displacement 

ber 

bol 

Increment 

Increment 

N 

T 

M 

A6n 

Afiip 

(kg) 

(kg) 

(cm-kg) 

(cm) 

(cm) 

401 

3zd 

19.8 

22.0 

3.9 

-0.00445 

0.01270 

402 

4zd 

35.8 

39.2 

7.3 

-0.00427 

0.01270 

403 

5zd 

51.8 

54.8 

23.5 

-0.00763 

0.02540 

404 

6zm 

67.8 

69.6 

46.4 

-0.00340 

0.01270 

405 

7zd 

83.8 

81.6 

45.4 

-0.00320 

0.01270 

421* 

4zm 

35.8 

37.8 

18.6 

-0.00358 

0.01270 

422* 

4zd 

35.8 

38.4 

17.4 

-0.00394 

0.01270 

423* 

4zd 

35.8 

37.3 

11.7 

-0.00374 

0.01270 

424* 

4zd 

35.8 

36 . 6 

20.2 

-0.00406 

0.01270 

425* 

6zd 

67.8 

69.3 

22.6 

-0.00774 

0.02540 

426* 

6zd 

67.8 

66.6 

54.1 

-0.00519 

0.02540 

427* 

5zd 

51.8 

51.4 

26.2 

-0.00722 

0.02540 

428* 

5zd 

51.8 

52.6 

20.6 

-0.00722 

0.02540 

429* 

7zd 

83.8 

81.6 

46.4 

-0.00765 

0.02540 

430* 

7zd 

83.8 

76.2 

14.3 

-0.00735 

0.02540 

431 

4zl 

35.8 

27.4 

9.7 

-0.00310 

0.02540 

432 

6zm 

67.8 

51.6 

19.8 

-0.00364 

0.02540 

433* 

4zl 

35.8 

26.4 

14.2 

-0.00119 

0.02540 

434* 

6zl 

67.8 

47.8 

23.8 

-0.00184 

0.02540 

435* 

6zd 

67.8 

66.5 

38.9 

-0.00770 

0.02540 

436* 

6zd 

67.8 

67.0 

35.8 

-0.00762 

0.02540 

437* 

6zd 

67.8 

68.0 

32.3 

-0.00750 

0.02540 

438* 

6zd 

67.8 

66.4 

30.7 

-0.00855 

0.02920 

439* 

6zd 

67.8 

65.4 

11.7 

-0.00712 

0.02540 

440* 

5zd 

51.8 

51.6 

18.5 

-0.00695 

0.02540 

441* 

7zd 

83.8 

79.9 

40.0 

-0.00717 

0.02540 

442 

3pd 

19.8 

20.2 

-7.3 

-0.00214 

0.00635 

443 

4pd 

35.8 

33.2 

-10.9 

-0.00339 

0.01270 

444 

5pd 

51.8 

44.5 

-31.4 

-0.00169 

0.00635 

445 

6pd 

67.8 

57.0 

-39. 3 

-0.00196 

0.00635 

446 

6nd 

67.8 

69.6 

66.0 

-0.00233 

0.00635 

447 

5nd 

51.8 

53.9 

53.0 

-0.00246 

0.00635 

448 

4nd 

35.8 

39.9 

33.7 

-0.00298 

0.00635 

449 

3nd 

19.8 

24.2 

17.7 

-0.00300 

0.00635 

450 

3zd 

19.8 

21.5 

6.6 

-0.00170 

0.00508 

451 

4zd 

35.8 

36 . 6 

13.1 

-0.00113 

0.00508 

Rotation 

Increment 

Ai|> 

(rad) 

-0.00113 

-0.00107 

-0.00243 

-0.00071 

-0.00088 

-0.00112 

-0.00092 

-0.00095 

-0.00082 

-0.00164 

-0.00123 

-0.00189 

-0.00187 

-0.00099 

-0.00175 

-0.00061 

-0.00073 

-0.00083 

-0.00049 

-0.00116 

-0.00125 

-0.00124 

-0.00174 

-0.00161 

-0.00162 

-0.00126 

-0.00053 

-0.00146 

-0.00055 

-0.00072 

-0.00003 

-0.00002 

-0.00007 

+0.00004 

-0.00022 

-0.00035 


*  tests  which  included  internal  cycling  or  non-typical  loading  path 


APPENDIX  III.  -  TABLE  III  (continued) 


DISPLACEMENT  INCREMENTS  AND  FORCES  AT  FAILURE 


Test 

Group  Normal 

Shear 

Moment 

Vertical 

Horizontal 

Num¬ 

Sym¬ 

Force 

Force 

Displacement  Displacement 

ber 

bol 

Increment 

Increment 

N 

T 

M 

(cm) 

A6t 

(kg) 

(kg) 

(cm- kg) 

(cm) 

452 

5zd 

51.8 

52.5 

15.3 

-0.00174 

0.00508 

453 

6zd 

67.8 

63.8 

22.2 

-0.00139 

0.00508 

454 

3nd 

19.8 

23.5 

18.2 

-0.00224 

0.00508 

455 

4nd 

35.8 

37.5 

33.4 

-0.00196 

0.00508 

456 

5nd 

51.8 

53.4 

47.6 

-0.00181 

0.00508 

457 

6nd 

67.8 

69.4 

69. 0 

-0.00186 

0.00508 

458 

3pd 

19.8 

20.3 

-7.9 

-0.00128 

0.00508 

459 

4pd 

35.8 

33.4 

-17.8 

-0.00152 

0.00508 

460 

5pd 

51.8 

43.5 

-26.4 

-0.00140 

0.00508 

461 

6pd 

67.8 

58.4 

-24.9 

-0.00131 

0.00508 

462 

lzd 

7.8 

9.4 

2.0 

-0.00195 

0.00508 

463 

2zd 

11.8 

13.1 

2.4 

-0.00190 

0.00508 

464 

3pm 

19.8 

14.4 

-4.1 

-0.00088 

0.00508 

465 

4pm 

35.8 

24.8 

-11.3 

-0.00058 

0.00508 

466 

5pl 

51.8 

32.6 

-33.6 

-0.00017 

0.00508 

467 

6pl 

67.8 

34.6 

-51.8 

+0.00035 

0.00508 

468 

3zl 

19.8 

13.1 

5.2 

+0.00008 

0.00508 

469 

4zl 

35.8 

24.9 

11.8 

-0.00010 

0.00508 

470 

3zl 

19.8 

15.0 

2.8 

-0.00366 

0.02540 

471 

4zl 

35.8 

25.1 

4.2 

-0.00198 

0.02540 

472 

5zl 

51.8 

37.5 

10.0 

-0.00099 

0.01270 

473 

6zl 

b7. 8 

45.6 

21.8 

-0.00047 

0.01270 

474 

3pm 

19.8 

15.5 

-13.4 

-0.00528 

0.02540 

475 

4pm 

35.8 

28.6 

-20.6 

-0.00493 

0.02540 

476 

5pl 

51.8 

32.2 

-24.4 

+0.00020 

0.01270 

477 

6pm 

67.8 

50.2 

-28.0 

-0.00470 

0.02540 

478 

3nm 

19.8 

16.9 

15.4 

-0.00226 

0.01270 

479 

4nd 

35.8 

32.2 

36.8 

-0.00694 

0.02540 

480 

5nl 

51.8 

39.6 

50.8 

-0.00180 

0.01270 

481 

6nm 

67.8 

57.5 

66.5 

-0.00222 

0.01270 

482 

2zm 

11.8 

11.5 

3.0 

-0.00288 

0.01270 

483 

3zm 

19.8 

17.0 

5.1 

-0.00206 

0.01270 

484 

4zm 

35.8 

30.3 

13.8 

-0.00226 

0.01270 

485 

5zm 

51.8 

43.5 

25.6 

-0.00187 

0.01270 

486 

6zm 

67.8 

56.6 

30.0 

-0.00189 

0.01270 

490 

6nd 

67.8 

64.7 

55.6 

-0.00503 

0.01750 

Rotation 

Increment 

A<jj 

(rad) 

-0.00029 

-0.00029 

-0.00003 

-0.00005 

-0.00005 

-0.00004 

-0.00038 

-0.00044 

-0.00045 

-0.00057 

-0.00041 

-0.00039 

-0.00036 

-0.00045 

-0.00038 

-0.00033 

-0.00024 

-0.00014 

-0.00105 

-0.00111 

-0.00033 

-0.00030 

-0.00204 

-0.00205 

-0.00111 

-0.00219 

+0.00017 

-0.00015 

+0.00013 

+0.00008 

-0.00092 

-0.00092 

-0.00088 

-0.00074 

-0.00063 

-0.00012 


APPENDIX  III.  -  TABLE  III  (concluded) 
DISPLACEMENT  INCREMENTS  AND  FORCES  AT  FAILURE 


Test 

Group  Normal 

Shear 

Moment 

Vertical 

Horizontal 

Rotation 

Num¬ 

Sym¬ 

Force 

Force 

Displacement 

Displacement 

Increment 

ber 

bol 

Increment 

Increment 

N 

T 

M 

(cm) 

A5 

A  <J> 

(kg) 

(kg) 

(cm-kg) 

(cm) 

(rad) 

491 

6zd 

67.8 

69.0 

27.6 

-0.00811 

0.02680 

-0.00135 

492 

6pd 

67.8 

66.1 

20.2 

-0.00562 

0.02540 

-0.00296 

493 

5nd 

51.8 

51.4 

■  37.3 

-0.00637 

0.02110 

-0.00052 

494 

5zm 

51.8 

51.1 

29.1 

-0.00575 

0.02670 

-0.00169 

495 

5pd 

51.8 

48.8 

-5.0 

-0.00398 

0.01640 

-0.00163 

496 

4zd 

35.8 

37.4 

11.3 

-0.00405 

0.01270 

-0.00041 

497 

4zm 

35.8 

37.6 

20.8 

-0.00296 

0.01270 

-0.00079 

498 

4pd 

35.8 

32.6 

-12.1 

-0.00140 

0.00635 

-0.00078 

499 

3nm 

19.8 

20.6 

14.4 

-0.00360 

0.01270 

-0.00053 

500 

3zd 

19.8 

21.6 

3.5 

-0.00374 

0.01270 

-0.00093 

501 

3pd 

19.8 

20.5 

0.4 

-0.00322 

0.01270 

-0.00152 

502 

6nm 

67.8 

51.0 

49.2 

-0.00179 

0.01270 

+0.00003 

503 

6zm 

67.8 

49.0 

5.8 

-0.00132 

0.01270 

-0.00034 

504 

6pl 

67.8 

47.4 

-12.5 

-0.00108 

0.01270 

-0.00092 

505 

5pl 

51.8 

39.1 

41.8 

-0.00145 

0.01270 

+0.00008 

506 

5zl 

51.8 

38.2 

23.0 

-0.00090 

0.01270 

-0.00030 

507 

5pl 

5l.8 

34.8 

-19.1 

-0.00060 

0.01270 

-0.00097 

508 

4nl 

35.8 

26.2 

26.8 

-0.00269 

0.02540 

+0.00016 

509 

4zl 

35.8 

26.3 

12.5 

-0.00118 

0.01270 

-0.00028 

510 

4pl 

35.8 

24.6 

-5.1 

-0.00017 

0.01270 

-0.00092 

511 

3nl 

19.8 

14.9 

13.7 

-0.00256 

0.02540 

+0.00009 

512 

3zl 

19.8 

14.9 

5.1 

-0.00254 

0.02540 

-0.00045 

513 

3pl 

19.8 

14.4 

-2.8 

-0.00214 

0.02540 

-0.00205 

514 

6zm 

67.8 

55.5 

16.3 

-0.00348 

0.02540 

-0.00192 

515 

5zm 

51.8 

43.5 

20.6 

-0.00178 

0.01270 

-0.00072 

516 

4zm 

35.8 

32.0 

9.1 

-0.00384 

0.02540 

-0.00151 

517 

3zm 

19.8 

17.5 

5.3 

-0.00440 

0.02540 

-0.00150 

518 

6nm 

67.8 

57.3 

34.6 

-0.00404 

0.02540 

-0.00075 

519 

5nm 

51.8 

46.0 

46.2 

-0.00571 

0.02540 

-0.00064 

520 

4nm 

35.8 

30.3 

28.2 

-0.00508 

0.02540 

-0.00051 

521 

3nm 

19.8 

18.0 

11.6 

-0.00232 

0.01270 

-0.00046 

522 

6pm 

67.8 

52.9 

14.0 

-0.00103 

0.01270 

-0.00103 

523 

5pm 

51.8 

43.3 

19.1 

-0.00176 

0.01270 

-0.00093 

524 

4pm 

35.8 

29.0 

-1.6 

-0.00164 

0.01270 

-0.00112 

525 

3pm 

19.8 

18.2 

-2.0 

-0.00211 

0.01270 

-0.00109 

APPENDIX  III.  -  TABLE  IV 


Means  and  Standard  Deviations  of 
Angles  Defined  by  Equation  (21) 


Angle 

Mean 

Value 
(deg. ) 

Standard 
Deviation 
(deg. ) 

0A 

39.02 

±3.66 

0B 

12.80 

±5.53 

ec 

26.76 

±3.17 

0D 

-3.09 

±1.88 

6E 

-4.32 

±3.33 

eF 

-1.23 

±2.64 

SG 

-38.79 

±3.80 

-12.07 

±6.05 

0T 

-26.72 

+3.21 

»*>»»** ’S’ 


Figure  1 
Figure  2 

Figure  3 
Figure  4 
Figure  5 

Figure  6 
Figure  7 

Figure  8 

Figure  9 

Figure  10 

Figure  11 

Figure  12 

Figure  13 

Figure  14 

Figure  15 


CAPTIONS  FOR  FIGURES 

Yield  Surface,  Stress  Path,  and  Plastic  Strain  Increment  Vector. 

Vertical  Cross-sectional  View  of  Direct  Shear  Device  in  the 
Plane  of  the  Applied  Tractions. 

Enlarged  View  of  a  Portion  of  Figure  2  Illustrating  Dimensions. 

Plan  View  of  Direct  Shear  Apparatus. 

Initial  Configuration  of  Soil  Specimen  and  Positive  Directions 
of  Applied  Tractions  and  Resulting  Displacements. 

Initial  &  Displaced  Configurations  of  Shear  Box  Device. 

Two  Dimensional  Projections  of  Loading  Paths  for  Tests  452  and 
472. 

Two  Dimensional  Projections  of  Displacement  Paths  for  Tests 
452  and  472. 

Graphs  of  Tangential  Force  Versus  Tangential  Displacement  and 
Normal  Force  Versus  Normal  Displacement  for  Test  452  and  472. 

Graph  of  Tangential  Force  Versus  Relative  Density  and  Current 
Tilt  Angle  for  Tests  452  and  472. 

Plots  of  Failure  Force  Coordinates  and  Displacement  Increment 
Vector  Projections  for  Tests  Having  Final  Normal  Loads  of 

19.8  kg.  and  35.8  kg. 

Plots  of  Failure  Force  Coordinates  and  Displacement  Increment 
Vector  Projections  for  Tests  Having  Final  Normal  Loads  of 

51.8  kg.  and  67.8  kg. 

Plot  of  Failure  Force  Coordinates  and  Displacement  Increment 

Vector  Projections  for  Tests  Having  a  Moment  at  Failure  between 
-10.1  cm. -kg.  and  -30.0  cm. -kg. 

Plot  of  Failure  Force  Coordinates  and  Displacement  Increment 

Vector  Projections  for  Tests  Having  a  Moment  at  Failure  between 
+9.9  cm. -kg.  and  -10.0  cm. -kg. 

Plot  of  Failure  Force  Coordinates  and  Displacement  Increment 

Vector  Projections  for  Tests  Having  a  Moment  at  Failure  between 
+29.9  cm. -kg.  and  +10.0  cm. -kg. 

Plot  of  Failure  Force  Coordinates  and  Displacement  Increment  Vector 
Projections  for  Test  Having  a  Moment  at  Failure  between  +49.9 
cm. -kg.  and  +30.0  cm. -kg. 


Figure  16 


CAPTIONS  FOR  FIGURES 


Figure  17  Pictoral  Identification  of  Vectors  and  Angles  Defined  by 

Equations  (19)  through  (21). 

Figure  18  Scatter  Diagram  of  the  Angle  0p  Versus  the  Relative  Density 

at  Failure. 

Figure  19  Scatter  Diagram  of  the  Angle  0p  Versus  the  Angle  of  Tilt  at 

Failure. 

Figure  20  Scatter  Diagram  of  the  Angles  0^  and  0 ^  Versus  the  Relative 

Density  at  Failure. 

Figure  21  Principal  Stress  Space  and  Failure  Theories. 


YIELD  SURFACE 


FIG.  1  YIELD  SURFACE,  STRESS  PATH,  AND  PLASTIC  STRAIN  INCREMENT 
VECTOR. 


FIG.  2  VERTICAL  CROSS-SECTIONAL  VIEW  OF  DIRECT  SHEAR  DEVICE  IN  THE 
PLANE  OF  THE  APPLIED  TRACTIONS. 


.  3  ENLARGED  VIEW  OF  A  PORTION  OF  FIG.  2  ILLUSTRATING  DIMENSIONS. 


PLAN  VIEW  OF  DIRECT  SHEAR  APPARATUS. 


FIG.  5  INITIAL  CONFIGURATION  OF  SOIL  SPECIMEN  AND  POSITIVE  DIRECTIONS 
OF  APPLIED  TRACTIONS  AND  RESULTING  DISPLACEMENTS. 


FIG.  6  INITIAL  £  DISPLACED  CONFIGURATIONS  OF  SHEAR  BOX  DEVICE. 
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FIG.  7  TWO  DIMENSIONAL  PROJECTIONS  OF  LOADING  PATHS  FOR  TESTS  452  AND 
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FIG.  8  TWO  DIMENSIONAL  PROJECTIONS  OF  DISPLACEMENT  PATHS  FOR  TESTS  452 
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FIG.  9  GRAPHS  OF  TANGENTIAL  FORCE  VERSUS  TANGENTIAL  DISPLACEMENT  AND 
NORMAL  FORCE  VERSUS  NORMAL  DISPLACEMENT  FOR  TEST  452  AND  472. 
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FIG.  12  PLOTS  OF  FAILURE  FORCE  COORDINATES  AND  DISPLACEMENT  INCREMENT 
VECTOR  PROJECTIONS  FOR  TESTS  HAVING  FINAL  NORMAL  LOADS  OF  51.8  KG.  AND 
67.8  KG. 


FIG.  13  PLOT  OF  FAILURE  COORDINATES  AND  DISPLACEMENT  INCREMENT  VECTOR 
PROJECTIONS  FOR  TESTS  HAVING  A  MOMENT  AT  FAILURE  BETWEEN  -10.1  CM. -KG. 
AND  -30.0  CM. -KG. 


FIG.  14  PLOT  OF  FAILURE  FORCE  COORDINATES  AND  DISPLACEMENT  INCREMENT 
VECTOR  PROJECTIONS  FOR  TESTS  HAVING  A  MOMENT  AT  FAILURE  BETWEEN  +9.9 
CM. -KG.  AND  -10.0  CM. -KG. 


FIG.  15  PLOT  OF  FAILURE  FORCE  COORDINATES  AND  DISPLACEMENT  INCREMENT 
VECTOR  PROJECTIONS  FOR  TESTS  HAVING  A  MOMENT  AT  FAILURE  BETWEEN  +29.9 
CM. -KG.  AND  +10.0  CM. -KG. 
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FIG.  16  PLOT  OF  FAILURE  FORCE  COORDINATES  AND  DISPLACEMENT  INCREMENT 
VECTOR  PROJECTIONS  FOR  TEST  HAVING  A  MOMENT  AT  FAILURE  BETWEEN  +49.9 
CM. -KG.  AND  +30.0  CM. -KG. 
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FIG.  18  SCATTER  DIAGRAM  OF  THE  ANGLE  0  VERSUS  THE  RELATIVE  DENSITY 
AT  FAILURE. 
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FIG.  21  PRINCIPAL  STRESS  SPACE  AND  FAILURE  THEORIES. 


